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ABSTRACT 


A  program  for  control  system  modeling,  simulation,  and 
realization  useable  by  the  novice  programmer  was  developed 
for  interactive  use  on  a  microcomputer.  The  program  uses  two 
software  packages;  PL/1-80  from  Elgltal  Research  and  EBASE  II 
from  Ashton  Tate. 

The  paper  covers  a  complete  example  for  modeling  and 

simulation  of  a  minimum  time  response  ripple  free  controller 
2 

for  a  1/S  Plant.  The  program  uses  two  sample  rates,  one  for 
modeling  the  computer  and  the  other  for  modeling  the  plant. 
Variable  delays  due  to  computer  computation  are  includea  and 
corrected  semi-transparent  to  the  control  designer. 

This  thesis  covers  two  different  fields  of  Interest,  one 
for  the  experienced  control  system  designer  yet  a  novice 
programmer,  and  the  other  for  the  novice  control  designer  yet 
an  experienced  computer  programmer.  Key  programming  concepts 
include  using  a  business  database  to  automatically  program  a 
control  problem  requiring  a  scientific  language. 
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I.  iNTsorysiigN 


A.  CONCEPT 

The  intent  of  this  thesis  was  to  develop  a  coir.puter 
program  that  an  experience!  controls  engineer  conid  use  to 
help  him  model,  simulate,  and  realize  control  systems  with 
the  convenience  of  a  desktop  microcomputer.  The  control 
engineer  will  not  need  to  he  an  experienced  computer 
programmer.  In  order  to  accomplish  the  above  requirements  it 
Is  necessary  that  the  program  be  easy  to  use,  effectively 
requiring  that  the  program  write  the  computer  code  Itself 
with  little  help  from  the  control  desiener.  This  is  commonly 
referred  to  as  automatic  programming.  In  order  for  t're 
program  to  be  a  useable  tool  for  the  controls  designer  the 
system  must  also  produce  a  good  model  of  the  s/stsm  tiat  is 
being  designed.  Obviously  all  control  systems  cannot  be 
considered  in  this  type  of  project.  The  approach  take'-  in 
this  thesis  was  to  use  control  systems  that  have  error 
signals  designed  to  be  equal  to  zero  or  differ  from  zero  by  a 
constant.  Other  guidelines  of  this  project  Include:  the 
equations  must  be  simple  difference  equations,  the  system 
must  allow  for  calculation  delays  within  the  computer  and 
these  delays  must  be  corrected  for  automatically  or  at  least 
allow  for  the  possibility  that  the  designer  ran  correct  for 
the  delays,  and  the  program  must  be  able  to  model  the  plant 


c 


between  the  sample  periods  of  the  computer  which  basically 
means  there  are  two  sample  periods,  one  for  the  computer  and 
one  for  the  plant.  Since  the  plant  cannot  be  programmed  as 
analog  it  must  be  discretized.  Uslnt,  a  faster  sample  rate  for 
the  plant  allows  simpler  equations  yet  maintains  a  good 
representation  for  the  plant,  iditionall.',  the  program  must 
be  able  to  model  analog  to  digital  and  lioital  to  analog 
conversions.  After  the  model  has  teen  designed  and  the 
simulation  run,  the  results  that  are  available  from  the 
program  must  be  easy  to  use  and  easy  to  interpret  by  the 
controls  engineer.  In  order  to  meet  all  of  these 
requirements,  most  of  the  program  development  itself  must  be 
transparent  to  the  controls  engineer.  The  approach  taken  was 
for  the  control  engineer  to  provide  only  equations,  function 
names,  and  function  relationships  by  provialng  the  names  -f 
other  functions  used  as  inputs  for  each  equation. 

B.  iSSO^'FTION3  AN!  RSCUIHE'1E:JT3 

The  assumptions  and  requirements  are  stated  below,  "irst, 
as  mentioned  above,  the  program  must  oe  relatively  easy  tc 
use.  The  program  must  have  high  piecision  math  capability 
since  it  was  designed  to  be  used  with  control  syst‘='*s.  Speed 
of  compiling  and  oneratioa  was  not  -onsiiered  criti'al  at 
this  point  since  microprocessors  are  becomlnr.  srallar  and 
faster  which  will  solve  any  realization  problems  i’-  the  near 
future.  A  major  consideration  is  that  once  *he  simulations 
have  been  run  and  the  results  are  acceptable  to  the  designer 
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then  the  simulation  equations  should  he  easily  converted  into 
realization  equations.  Basically  this  means  the  simulation 
equations  in  the  computer  should  he  easily  converted  into 
hardware  equations  that  can  he  used  to  control  a  real  plant. 
In  order  to  accomplish  this  task  it  was  decided  that  the 
program  must  develop  source  code  that  could  he  compiled  ty  a 
high  level  language  compiler  for  either  computer  simulation 
or  hardware  implementation.  This  approach  makes  the  prrf'’'am 
transportable  and  hardware  independent. 

C.  IMPLEMENTATION 

This  project  was  implemented  entirely  on  microprocessor 
based  systems  to  show  that  it  is  feasiole  and  practical  to 
use  these  small  systems  for  major  development.  The  present 
system  is  slightly  constrained  due  to  the  length  cf  compile 
time.  All  programs  were  developed  using  a  Z8ic  microprocessor; 
nowever,  the  code  is  written  in  high  level  source  languages 
therefore  the  programs  will  also  run  on  the  simpler  S?Si 
based  microcomputers  or  on  the  more  complex  sixteen  bit 
microprocessors  such  as  the  6066. 

The  languages  chosen  were  PL/l-60  from  ligital  Research 
oecause  of  the  math  and  hardware  capabilities  f  the  langua'?*^ 
allowing  for  control  of  analog  to  digital  ani  other  nrooesses 
inside  the  computer,  and  DBASE  II  from  Ashton  Tate  .hi'h 
allows  easy  and  high  level  -manipulation  of  -^.atabases  and  file 
systems.  Cne  cf  the  major  points  of  this  thesis  is  to 
manipulate  functions  and  functional  i'elation<ihins  as 
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databases.  Once  the  relationships  have  been  developed,  the 
database  program  writes  PL/1  source  cole  from  the  database 
information.  The  PL/1  source  code  is  then  compiled  by  the 
PL/l  compiler  and  the  simulation  run  on  the  resident 
hardware. 

The  present  programs  and  the  test  example  In  this  report 
provide  simulation  only  and  send  the  results  in  floating 
point  format  to  a  disk  file.  The  floating  point  file  can  then 
be  converted  to  display  the  information  in  tabular  or 
graphical  form  such  as  shown  at  the  back  of  this  report.  The 
information  can  be  displayed  in  a  variety  of  forms  such  as  on 
display  consoles,  printers,  plotters,  or  graphic  terminals. 
The  real  power  for  design  is  in  the  use  of  graphical 
equipment  to  provide  the  control  engineer  with  the  entire 
picture  of  the  system  under  test.  The  Hewlet  Packard  9872B 
plotter  and  2647A  graphics  terminal  and  the  Intercolor  8363 
graphics  terminal  were  used  to  show  that  the  data  produced 
was  basically  machine  independent.  The  cole  can  ’le  compiled, 
run,  and  displayed  on  a  variety  of  different  machines. 


II.  GENERAL  EISCUS3I0N 


A.  CONTROL  ENGINEER  INFORMATION 

The  control  system  engineer  does  not  need  to  know  how  the 
computer  programs  work  In  order  to  use  the  system.  It  is 
assumed  that  the  control  engineer  is  a  novice  programmer  with 
some  knowledge  of  the  CP/M  operating  system  and  limited 
background  in  program  writing;  however,  it  is  assumed  he  is  a 
Knowledgeable  control  designer.  This  means  be  is  fairly  well 
versed  in  taking  control  systems  and  producing  discrete 
system  equations  from  the  continuous  control  equations.  These 
equations  are  designed  and  optimized  by  whatever  means  a-e 
available  either  by  hand  or  on  a  mainframe  computer. 
(Microcomputers  will  soon  be  able  to  do  major  system 
development  equations.)  Once  the  equations  are  developed  the 
designer  enters  them  into  a  microcomputer  and  tests  the 
simulations.  If  the  simulations  are  acceptable  then  it  should 
be  an  easy  matter  to  convert  the  simulation  code  for  the 
controller  design  into  development  cede  for  a-'t'cal 
implementation  into  hardware. 

For  ease  of  use  the  programs  are  menu  driven  and  use 
simple  cursor  arrows  or  commanus  for  pointing  at  selections. 
Additionally  very  basic  information  is  needed  such  as  the 
name  of  functions  and  equations  in  a  simple  format  the 
computer  can  use.  This  procedure  is  accomplished  oy  using  the 
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enter 


database  program  which  only  allows  a  user  to 
information  In  very  specific  fields  of  very  specific  type. 
The  type  of  information  required  from  the  designer  is  to  add 
or  delete  equations,  to  change  initial  conditions,  and  to 
provide  function  names. 

There  is  an  initial  learning  curve  with  this  program  but 
an  attempt  has  been  made  to  keep  it  as  short  as  possiole. 
There  are  some  peculiarities  put  in  as  a  requirement  to  make 
interfacing  between  the  human  and  the  computer  posslbl«>.  Some 
examples  include  names  such  as  INPUT :A  or  INPUT:E  that 
indicate  which  function  the  designer  is  talking  accut  to  the 
computer.  As  an  example  INPUT:A  will  associate  a  function  in 
the  computer  with  the  function  name  the  designer  has  selected 
such  as  INPUT:A  *  Un.  Another  very  abstract  name  is  AIPPP. 
The  user  Is  expected  to  enter  either  1  if  he  wants  the 
PRESENT  value  of  the  function  named  under  INPUT  :A  or  2  if  he 
wants  the  EELATEE  value  of  the  function  named  under  INPUT  :A. 
If  the  designer  will  be  using  both  the  present  and  past 
values  then  enter  1  to  get  the  present  value  and  the  oast 
value  will  also  be  available.  The  strange  requlremert  for 
A1P2D  Is  provided  only  to  allow  much  faster  programs  when  th“ 
code  is  implemented  in  hardware  and  the  designer  only  needs 
past  or  delayed  values  of  a  function.  Two  other  codes  are  lAP 
and  lAE  for  INPUT  A  PRESENT  and  INPUT  A  EELAYEE.  These  two 
codes  are  provided  to  keep  the  equations  needed  for  the 
functions  short.  Another  variable  is  "o"  which  is  asec  tc 


define  the  delayed  or  previous  time  sample  of  the  functions 
own  OUTPUT.  There  are  similar  definitions  for  up  to  five 
inputs  INPUT:A  through  INPUT:E  for  each  function,  k  function 
can  he  defined  in  three  ways;  "F"  means  the  function  is  an 
equation  that  the  computer  finds  a  register  to  bold  the 
result  in,  '*ir(3c)"  means  the  value  is  to  be  placed  in  one  of 
eight  analog  to  digital  registers,  and  "ikix)"  means  the 
value  is  to  be  placed  In  one  of  eight  digital  to  analog 
reclsters.  An  example  equation  and  simple  block  ilagrafr  is 
shown  in  Figure  1. 

Other  variables  used  Include  *t"  for  the  discrete  sample 
period  and  "Ti"  for  the  sample  period  of  the  plant.  If  the 
value  of  T  is  normalized  and  Tl  has  the  same  value  as  T  then 
the  programs  will  run  faster  in  the  realization  of  the  code 
since  these  values  effectively  disappear  from  the  equations. 
It  should  be  noted  that  all  functions  are  done  in  floating 
point  so  the  simpler  the  equations  the  faster  the  code.  The 
example  problem  at  the  end  of  this  paper  normalizes  T  but 
still  includes  it  in  the  equations  to  show  how  it  could  show 
up  in  the  function.  Finally  the  designer  is  allowed  to 
specify  end  criteria  and  intervals  for  calculating  the  plant 
relative  to  the  computer  Itself, 

After  the  designer  has  entered  the  equations  into  the 
database  and  is  satisfied  that  all  function  relationships  are 
correct  then  all  that  is  required  is  to  enter  the  correct 
menu  selection  that  compiles  the  code.  The  ccmuuter  program 
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Assume  by  previous  deflultlon  that  INPUT:A  «  X.  INPUT  :B  »  T 
FUNCNAME  «  Uz ,  and  the  below  equatioo  is  desired: 

D(k)  -  X(k)  -  (2  ♦  X(k-l))  ♦  y(k)  -  (0.b  *  D(k-l)) 

Making  substitutions: 

F  »  Ez  *  E(k),  lAF  »  X  »  X(k),  lAD  »  X(k-l),  IBP  »  I  »  T(k), 

and  0  *  OUTPUT  EELATED  »  D(k-l)  leads  to  the  equation  - 

F  »  lAP  -  (2  *  lAD)  ♦  IBP  -  (a. 5  *  0) 

Notice  that  the  above  function  only  has  two  inputs,  X  ann  T 


Figure  1 

Simple  Difference  Equation  leTionstrating  System  Variables 


will  oov  write  the  necessary  source  code  from  the  database 
and  subsequently  invoice  the  PL/1  compiler  which  will  turn  the 
code  into  a  program  that  will  run  on  any  of  the  previously 
discussed  CP/M  based  systems.  The  executable  file  produced  is 
called  ANSWER. COM  and  is  run  on  a  CP/M  based  system  by  typing 
"answer"  in  response  to  the  operating  system  input  prompt. 
Luring  program  ANSWER. COM  execution  the  results  will  be  sent 
to  a  floating  point  file  called  LATA.PLT.  Also,  while  the 
program  is  running  the  Intermediate  results  are  printed  on 
the  display  device  for  program  troubleshooting  should  the 
results  not  be  as  expected.  The  file  EATA.FLT  is  not 
printable;  however,  an  example  utility  program  called 
PRINT. PLI  is  provided  that  changes  the  floating  point  file 
into  a  printable  text  file.  This  program  is  executed  by 
typing  PRINT  PILE  where  FILE  is  the  name  of  the  designer 
provided  text  file  the  designer  wants  the  results  placed 
into.  With  further  programming  development  the  floating  point 
data  can  be  converted  into  any  format  needed  by  different 
aevices  such  as  plotters  as  shown  by  the  plot  results  in 
Appendix  C. 

The  information  for  understanding  the  database  program 
operation  is  provided  in  Chapter  V.  This  section  is  not 
necessary  for  understanding  at  the  design  level  but  it  does 
explain  variables  that  will  be  needed  in  equation 
development.  Additionally  the  example  problem  will  be  helpful 
to  the  design  engineer,  when  the  designer  is  ready  to  use  the 


program  all  that  is  needed  Is  to  enter  "CBASE  MENU"  to  the 
CP/M  operating  system  prompt.  As  previously  stated  the 
designer  will  need  some  familiarity  with  CP/M  before  the 
program  can  be  of  use. 

B.  PROGRAMMER  INFORMATION 

This  section  covers  the  experience  necessary  for  program 
maintenance  and  further  development.  These  programs  would  be 
worked  on  and  coded  by  an  individual  that  has  a  general 
control  system  background;  however,  he  must  have  substantial 
experience  with  computer  programming.  The  type  of  computer 
experience  needed  pertains  to  operating  systems ,  event 
counters,  and  automatic  programming.  The  programmer  must  have 
a  good  knowledge  of  both  PL/l-80  and  IBASE  II.  Additionally 
it  would  be  necessary  to  know  microcomputer  hardware 
implementations  such  as  analog  to  digital  converters  and 
sample  and  holds.  These  devices  are  simulated  in  the  computer 
programs  but  further  useful  extensions  would  involve 
implementing  the  code  into  realized  hardware.  The  information 
required  for  the  programmer  to  understand  the  programs  Is 
contained  in  the  remainder  of  this  thesis.  The  majority  of 
this  thesis  is  the  computer  code  required  to  implement  a 
design  environment  for  a  novice  programmer  through  the  use  of 
automatic  programming. 
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III.  SYSTEM  I£I0SXfiCR4SIES 


Major  problems  encountered  in  this  thesis  had  to  do  with 
programming  and  efficiency.  The  major  problem  Is  that  there 
is  no  one  programming  language  that  contains  all  the  features 
necessary  to  implement  the  entire  project.  LBASZ  II  was 
chosen  because  of  the  relational  database  characteristics 
even  though  It  has  very  limited  math  or  hardware 
capabilities.  The  requirements  for  reasonable  math  is  the 
reason  PL/l-30  was  chosen  and  though  some  database  systems 
are  available  for  PL/1  they  do  not  have  all  the  features  of 
IBiSS  II.  K  major  deficiency  of  XBASE  II  is  that  it  will  not 
allow  the  printing  of  a  semicolon.  The  semicolon  is 
considered  a  command  character  in  XBASE  II  and  therefore 
reservea.  This  can  be  considered  a  major  bug  for  a  database 
system,  especially  one  that  was  designed  for  business 
applications.  It  Is  foolish  to  oelieve  business  will  not  want 
to  include  a  semicolon  someplace  in  one  of  their 
transactions.  Ashton  Tate  was  notified  and  they  replied  it 
was  not  worth  trying  to  fix.  This  presents  a  problem  since 
the  database  program  was  used  to  write  the  source  code  in 
PL/1  and  PL/1  requires  a  semicolon  at  the  end  of  every 
statement.  The  problem  was  wortced  around  by  having  XBASE  II 
write  a  backslash  every  place  a  semicclon  was  needed  in  tae 
PL/1  source  text.  Next  c-  PL/1  preprocessor  called  SEFCR'^.PLI 


was  used  to  change  all  backslashes  to  semicolons  prior  to 
invoking  the  PL/1  compiler.  This  takes  a  large  amount  of  the 
compile  time  which  could  have  been  easily  corrected  by  a 
reasonable  response  from  Ashton  Tate.  PL/1  also  has  some 
problems;  however,  Eigital  Research  was  more  than  happy  to 
try  to  overcome  deficiencies.  There  are  several  problems  in 
their  floating  point  routines  due  mainly  to  conversion  rules 
which  can  cause  strange  effects  in  the  results.  Most  floating 
point  functions  such  as  the  exponential  function  tend  tc 
convert  the  arguments  to  fixed  before  calculation  and  this  in 
effect  truncates  the  arguments.  As  an  example  if  one  desires 
to  raise  an  exponent  to  a  power  and  that  power  is  a  function 
such  as  X/Y  then  both  X  and  Y  and  the  result  get  truncated. 
If  X/Y  is  a  fraction  less  than  1.0  the  effect  is  to  truncate 
the  result  to  e.0  and  the  result  of  the  exponent  to  a 
fractional  power  always  ends  up  equal  to  1.0.  The  way  around 
this  problem  is  to  convert  X  and  Y  to  CHAfiACTKRS  before  doing 
the  division.  This  Is  only  part  (f  the  problem,  slven  If  X  and 
Y  origionally  were  floating  print  and  were  converteo  to 
characters  the  result  of  the  exponent  evaluation  Is  *’lxed 
decimal.  Now  this  result  must  je  converted  by  another 
function  called  JTC  In  order  to  get  the  result  back  to  tr.e 
system  in  floating  point.  An  example  of  this  shows  up  In  the 
MONITOR. PLI  program  for  the  function  called  OF,  This  does 
cause  some  limitation  on  the  type  of  equations  that  the 
designer  can  use  without  knowing  some  peculiarities  of  the 
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designer  can  use  vlthout  knoving  sone  peculiarities  of  the 
system. 

For  the  sake  of  program  efficiency,  as  far  as  development 
feoes,  all  variables  or  functions  the  designer  uses  are 
considered  floating  point.  This  makes  programming  easier  with 
higher  precision;  however,  the  code  takes  longer  to  execute 
even  though  the  execution  time  was  not  considered  critical  in 
this  Implementation.  Using  floating  point  for  all  functions 
Introduces  some  Inaccurate  realizations  when  dealing  with 
analog  to  digital  and  digital  to  analog  converters.  This  is 
not  a  major  problem  since  these  devices  will  have  greater 
accuracy  in  the  future  or  the  devices  can  be  modeled  (tiore 
accurately  In  future  designs. 

All  results  presently  produced  by  this  system  are  placed 
on  floppy  disks  inside  the  computer.  Writing  to  the  iisK*  and 
displaying  the  results  on  a  terminal  are  an  obvious  slowdown 
in  system  operation.  These  Inefficiencies  do  not  show  up  in  a 
realization  of  the  control  system  since  the  controlle"  will 
be  self  contained  computer  code  communicating  only  with  the 
converters.  Additionally  the  code  will  be  faster  since  it 
will  not  be  necessary  to  simulate  the  plant. 

Presently  the  realization  of  the  cede  is  not  imple'^ented . 
The  software  code  needed  for  converters  and  the  actual 
monitor  needed  for  the  realized  hardware  is  not  ieveloT.ed. 
This  is  not  a  major  problem  and  should  be  relatively  straight 
forward  for  a  computer  programmer  with  general  haraware 
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experience.  The  actual  control  program  will  he  identical  and 
will  only  need  to  link  into  the  new  converter  programs.  In 
addition  the  simulation  display  is  not  fully  implemented.  To 
make  this  project  a  powerful  tool  the  simulation  results 
should  he  presented  in  graphical  form  to  the  designer  similar 
to  the  plots  at  the  end  of  this  taper.  If  properly  displayed, 
when  the  designer  changes  one  part  of  the  system  all  effects 
throughout  the  system  should  he  dynamically  evident.  The 
advantages  of  a  completely  contained  hardware,  software,  and 
display  system  are  obvious. 

The  time  required  to  compile  the  programs  is  definitely  a 
function  of  the  computer  hardware.  Most  newer  ZSZ  based 
single  user  systems  can  compile  the  code  in  about  fifteen 
minutes.  8080  based  systems,  older  260  based  systems,  and 
multiuser  systems  that  share  the  central  processor  will 
require  up  to  thirty  minutes  to  compile.  Recall  tae  compile 
time  is  not  critical  as  mentioned  earlier  ann  the  realization 
code  will  run  in  the  order  of  mllllsecoads .  As  a  final  note 
the  programs  would  run  much  faster  if  only  one  language  and 
no  preprocessor  were  needed. 

Another  area  that  requires  discussion  is  the  theory 
involved  in  the  design  of  the  correction  factor.  levelopnent 
of  a  general  prediction  factor  due  to  delays  caused  by 
conversions  and  computer  calculations  is  covered  very  briefly 
in  this  paper.  The  practical  results  of  a  predictor  that  is 
used  with  the  example  problem  is  shown  in  the  graphs  at  the 


E2 


IV.  FUTORi:  EXTSNSIONS 


After  a  careful  Investigation  of  the  types  of  protlerrs 
encountered  with  various  hardware  and  software  interfaces,  it 
becomes  very  evident  that  a  sixteen  hit  or  larger 
microprocessor  that  offers  high  speed  calculation  integrated 
into  a  complete  graphic  and  plotting  machine  will  he  a 
definite  boost  to  the  realization  of  future  work. 
Additionally  a  math  processor  integrated  circuit  that  can 
operate  concurrently  with  the  central  processor  will  be  of 
invaluable  Importance.  Once  you  place  these  hardware  items 
along  with  good  converter  circuits  into  one  unit  it  will  r.o 
longer  be  necessary  to  simulate  in  one  nachine  and  realize  in 
another.  Once  the  simulation  code  works  all  that  will  be 
needed  will  be  to  specify  executing  the  realization  code  on 
the  same  machine  if  connectors  are  available  on  the  machine 
for  communicating  with  the  plant.  Once  the  systerr  ooerates 
correctly  the  hardware  would  be  reduced  to  a  smaller  special 
purpose  system;  however,  the  power  of  using  a  microprocessor 
for  the  controller  cannot  be  overlooked.  If  one  provides  the 
capability  of  a  full  blown  microcomputer  in  the  zontrol  loop 
even  though  all  processes  are  not  needed,  the  system  could  be 
sent  to  a  distant  changing  environment  and  if  the  control 
equations  needed  to  be  changed  due  to  unforseen  ciroumstancfs 
it  could  be  done  remotely.  £y  using  the  above  method  it  would 
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be  possible  to  reconfigure  and  realize  a  new  system 
dynamically  which  cannot  be  done  by  a  physical  one  time 
hardware  implementation. 

mother  requirement  is  the  software  system  needed  for 
this  type  of  development.  The  problems  of  integrating  two 
completely  different  software  packages  is  evident  in  this 
paper.  4n  ideal  software  package  would  be  a  scientific 
relational  database  that  has  the  capaoility  to  control  the 
computer  hardware  and  provide  sell  modifying  code  or 
automatic  programming  for  it's  own  operating  system.  This 
type  of  language  would  make  the  programming  more  efficient, 
ease  the  implementation,  and  provide  a  friendlier  environment 
for  future  expansion.  This  type  of  software  package  is 
realizable  under  the  new  concurrent  languages  being  developed 
such  as  AEA  (product  of  the  tepartment  of  Eefence).  AIA  has 
not  been  proven  to  work  in  microcomputers  since  it  is  such  an 
extensive  language?  however,  several  companies  are  showing  it 
to  be  feasible  in  the  newer  computers  such  as  the  €84300 
microprocessor  from  Motorola.  The  66000  is  a  thirty  tw  :  bit 
processor  with  a  sixteen  bit  external  ous .  Motorola  is 
presently  designing  a  version  of  the  66000  tc  nave  a  thirty 
two  bit  external  bus. 

It  is  recommended  that  future  systems  contain  processors 
such  as  the  68000  and  math  processors  such  as  Intel's  806?  or 
similar  math  processors  from  other  companies  such  as  Tetas 
Instruments.  Additionally  high  speed  high  precision  analof,  tn 


digital  and  digital  to  analog  converters  should  be  used.  The 
type  of  converters  used  should  contain  their  own  multiple 
input  sample  and  hold  circuits.  This  would  insure  capturing 
all  inputs  at  the  same  time  (versus  time  multiplexing)  making 
the  control  realizations  more  practical  for  real  time 
operations. 


V.  EgSCRIPTION  0?  database  ALSQglTHM 

A.  (JENERAL  DISCUSSION  (Refer  to  Appendix  A) 

This  section  covers  how  the  database  turns  the  functions 
into  PL/1-80  code.  It  covers  each  suorcutine  and  there  ar“ 
hints  on  how  the  control  designer  uses  txie  systerr;.  The 
easiest  way  for  the  control  designer  to  use  this  thesis  is  to 
briefly  read  this  material,  go  throUf.h  the  example  pro  lem. 
and  then  proceed  directly  into  the  program.  In  general  a 
database  is  a  program  that  contains  related  information. 
Normally  these  urograms  are  used  in  business  applications 
where  it  is  desired  to  keep  track  of  an  indiviiuals  name, 
address,  or  account  information;  however,  databases  can  he 
used  to  store  other  information.  In  this  program 
implementation  the  database  contains  the  function  name,  the 
type  of  function  such  as  an  analog  function,  names  of  other 
functions  that  are  used  by  this  function  as  inputs.  the 
equations  for  the  functions,  the  relationship  cetween 
functions,  and  if  the  function  will  be  saved  for  later  r9>-all 
and  plotting. 

The  program  has  the  ability  to  develop  new  or  recall  old 
system  models  that  have  been  previously  developed.  Once  the 
program  is  entered  the  model  can  te  changed  by  aiding  or 
deletlnt,  transfer  functions,  change  relationsnips .  or  change 
the  equations.  Additionally  the  models  initial  conditions  cm 


be  changed.  If  only  the  initial  conditions  are  changed  the 
program  will  compile  much  faster  since  each  function  will  not 
need  to  be  recompiled.  After  the  system  model  is  developed 
the  database  proceeds  to  write  PL/l  "include"  files  that  are 
later  included  into  standard  PL/l  modules  developed  for  this 
project . 

The  features  of  the  LEASE  II  program  are  the  high  level 
relational  database  language  and  the  ability  of  the  program 
to  write  text  files  onto  floppy  disk.  The  disadvantages 
include  a  limited  math  capability  and  it  does  not  provide  the 
ability  to  interface  with  the  microcomputer  hardware.  A 
severe  drawback  in  the  program  is  that  it  cannot  produce  a 

semicolon  in  any  of  the  text  files  tnat  it  writes, 

B.  LEVEL  #1  (MAJOR  SOBHOUTINSS  AND  PROCEDURES} 

1*  MEmiiiCME  -  initial  menu  selection 

a.  Initializes  all  variables  and  the  database. 

b.  Enters  a  CASE  statement  to  select  a  new  problem, 

old  problem,  save  the  present  problem,  display  present 

function  names,  change  only  the  initial  conditions  of  the 

present  problem,  or  quit  the  program. 

c.  See  Figures  15A  and  15B  in  Appendix  A. 

2.  CMENU^CMJ  -  second  level  menu  selection 

a.  Develops  the  control  model  database. 

b.  Enters  a  CASE  statement  to  edit,  add,  delete, 
recall  deleted  functions,  or  change  initial  conditions. 


c.  Calls  CONTROL. CMI. 

d.  See  Figures  16A  and  16B  in  Appendix  A. 

3.  CONTROL^CMI  -  database  control  program 

a.  Calls  procedures  LCLMON.CWr,  PHOC.CME,  DATA.:*"!, 
and  FORM. CM!  that  develop  the  FL/1  "include”  source  code. 

b.  Ends  the  database  program  and  invokes  the  PL/1 
programs . 

"  monitor  declaration  program 

a.  develops  ECLKON.TXT  for  the  i^OMTOH.PLI  program. 

b.  Calls  MCNICL.CMr. 

2.  MONECL^CML  -  develops  AA.TXT,  AE.TXT,  ID. TXT,  ani 
EA.TXT  include  functions  for  '^CMTOR.PLI. 

6*  ”  Eevelops  ICLFUN.XXX,  INPUT. XXI,  FUNC.XXX, 

and  OUT. XXX  for  each  function  developed  by  the  control 
engineer  to  be  used  is  the  appropriate  function  AA.PLI, 
AD.PLI,  DD.PLI,  or  lA.PII. 

-  data  declaration  ?rocedu’'£ 

a.  Develops  EATAINFO.TXT  to  be  used  in  PRINT. ?LI 

b.  lATAlNFC  .TXT  contains  valuable  format  information 
pertaining  to  the  EATA.JLT  file. 

c.  Calls  DATAI.CMI  -  additional  data  format  file 

"  PL/1  command  erecutlor.  program 

a.  Eevelops  submit  file  PORT. SUB  used  to  control  the 
order  of  PL/1  program  compiling. 

b.  Issues  commands  to  reformat  files  containing  a 
backslash  Into  files  with  a  semicolon. 
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c.  Calls  CCLFUN.C^!£  function  name  setup  for  REFORM' 

d.  Calls  REFORM.CMD  -  general  format  for  PL/1 
compiling. 

9.  SHORTj.CMD  -  short  compile  program 

a.  levelops  initial  conditions  for  LCLKON.TXT  and 
MONITOR. PLI. 

b.  Issues  commands  to  control  the  compiling  of  only 
MONITOR. PLI  by  developing  a  submit  file  SfORM.CMI. 

c.  Calls  SFORM.CME  -  short  format  for  PL/1  compiling 
similar  to  FORM. SUB  called  SFORM.SUB. 

C.  LEVEL  #2  (PROGRAM  DESCRIPTION) 

Since  the  majority  of  the  code  is  in  higher  level 
languages,  and  is  designed  in  modules,  it  is  already  ir 
algorithmic  form.  It  is  assumed  that  the  individual  needing 
to  understand  the  code  will  be  familiar  with  references  1 
through  4;  therefore,  an  extensive  discussion  of  the  computer 
code  will  not  be  provided.  The  example  problem  will  cover  all 
the  features  of  the  program. 

The  basic  structure  of  the  database  Includes  a  memory 
file  that  contains  the  value  of  initial  variables  and  a 
database  file  that  contains  the  functions  and  their 
relationships  to  each  other.  The  entire  database  program  is 
involved  with  manipulating  and  recalling  information  from  the 
two  above  mentioned  files. 


The  initial  memory  file  is  shown  in  Figure  2  and  is 
contained  in  file  INIT.MEM.  This  file  is  provided  so  that 


initial  variables  called  will  not  produce  errors  in  the 
database  program.  The  main  variables  of  concern  are: 

1.  T  ~  discrete  sample  period, 

2.  T1  -  analog  sample  period, 

3*  PtiUil  ”  criteria  for  tne  digital  simulation, 

4.  ALIMIT  -  number  of  analog  samples  for  each  digital 
sample. 


5.  EELAY  -  delays  due  to  calculation  and  conversion 
expressed  as  the  number  of  analog  counts  after  the  sample 
period  T. 

Unseen  system  variables  include: 

1.  ACTION  -  users  next  menu  selection, 

2*  PARAM  -  holds  database  program  name, 

3.  CN  -  test  variable, 

4.  FN  -  test  variable, 

5.  TEMPOUT  -  temporary  variaole, 

6.  elements  -  data  array  position  for  a  variable, 
i!i£  "  index  counter, 

8.  FUN  -  string  of  characters  that  is  used  to  uniquely 
define  a  function  flletype, 

9*  LATAI  -  number  of  data  array  elements  saved. 

The  basic  structure  for  all  functions  is  shown  in  Figure 
3  and  is  contained  in  file  INIT.EBF.  The  initial  values  of 
the  database  are  shown  at  the  bottom  of  the  Figure  3.  The 
only  items  the  control  designer  sees  are  listed  under  the 
NAME  column  of  Figure  3.  The  user  must  enter  values  for  the 
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ACTION 

(C) 

0 

PAHAM 

(C) 

INIT 

CN 

(C) 

INIT 

FN 

(C) 

FUNCTION 

T 

(N) 

1 .000 

ELIMIT 

(N> 

10 

ALIMIT 

(N) 

5 

TSrPOOT 

(C) 

000 .0000 

ELEMEMTS 

(N) 

0 

IND 

(N) 

0 

FUN 

(C) 

AiClE?GHIJKLt^NOP0fi3TiJV«XYZ 

EELAY 

(n) 

0 

T1 

(N) 

1  .000 

EATAl 

(N) 

0 

Figure  2 

Variables  Used  for  Initialization 
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STRUCTURE  FOR  FILE:  INIT.EBF 
NUMBER  OF  RECORCS :  00001 

DATE  OF  LAST  UPDATE:  00/00/00 


PRIMARY 

USE  EAT  ABASE 

FLE 

NAME 

TYPE 

width 

001 

FUNCNAME 

C 

008 

002 

FUNCTYPB 

C 

002 

003 

NR: INPUTS 

C 

001 

004 

INPUT: A 

C 

008 

005 

AlP2t 

C 

001 

006 

INPUT  :B 

C 

008 

007 

B1P2L 

C 

001 

008 

INPUT.'C 

c 

008 

009 

C1P2D 

c 

001 

010 

INPUT:! 

c 

008 

011 

D1P2L 

c 

001 

012 

INPUT :E 

c 

008 

013 

E1P2E 

c 

001 

014 

FK 

N 

008 

015 

FEL 

N 

003 

016 

EQUATION 

c 

050 

017 

SAVEEATA 

L 

001 

018 

EATANAME 

C 

005 

*  Below  Is  a  list  of  the  values  used  in  the  initialization. 
These  values  correspond  directly  to  the  structure  above. 


00001 


new  func 

kl 

0 

ONE 

1 

TWO 

1 

THREE 

1 

FOUR 

1 

FIVE 

1 

0.0000 
0 .0000 
Ar(i)=0 

AD(1) 


Figure  3 

Initial  Eatabase  Values 


variables  under  NAME  of  the  correct  type  and  size  specified. 
The  definitions  for  the  variables  are  given  below: 

1.  FUNCNAMI  -  function  name  used  to  Identify  this 
function  initialized  to  "NEW  FUNC". 

2.  FDNCTYPE  -  describes  this  function's  type  as; 

a.  AA  -  analog, 

b.  AD  -  analog  to  digital, 

c.  EE  -  digital, 

d.  DA  -  digital  to  analog. 

3*  "  This  is  the  number  of  other  functions 

this  function  uses  for  inputs.  AL(x)  and  EA(x)  registers  are 
not  counted  as  functions. 

4.  INPUT^A  -  This  is  a  variable  for  name  of  function 
used  as  an  input  and  is  shown  at  tne  bottom  of  Figure  3  as 
"one".  This  value  is  referred  to  in  the  equations  as  lAP  or 
lAL  for  INPUT  A  PRESENT  or  INPUT  A  DLLAYEE. 

5*  A1P2D  -  This  is  a  code  to  tell  tne  computer  the 
designer  is  talking  about  the  "a'  input  and  he  wants  either  a 
"l"  for  the  PRESENT  value  or  "2"  for  the  lELiYEE  value. 
Entering  a  "l'  will  make  both  the  present  and  delayed  values 
available  but  the  program  will  run  slower  than  if  only  the 
aelayed  value  is  needed. 

6.  lN£UTiB  -  This  variable  is  Identical  to  4  aLcve 
except  this  is  the  "B”  input  and  has  an  initial  name  of  T.VO" 
as  shown  at  the  bottom  of  Figure  3.  INPUTiC,  IN?UT:E,  and 
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INPnT:E  are  the  same  as  above  which  allows  for  up  to  five 
inputs . 

'?•  11IS5  ”  This  variable  Is  the  same  as  5  above  except 

it  is  for  the  "B*  input  function.  C1P2E,  D1P2L,  and  E1P2I;  all 
have  the  same  function  except  they  are  used  with  their 
respective  Inputs. 

8.  -  This  function  is  the  present  value  of  the 

function  named  under  FUNCNAME  above. 

S*  “  this  is  the  delayed  or  previous  sample  period 

value  of  the  function  named  in  FUKCNAME  '.le.  FK  delayed). 

10.  EQUATION  -  this  is  the  equation  for  the  function 
named  in  fUNCNAME  above.  The  equation  takes  three  forms*' 

F  =  — — - this  is  a  function, 

AI(x)  *  -  value  for  an  analog  to  digital, 

DA(x)  - value  for  a  digital  t('  analog. 

11.  SAVEIATA  -  logical  value  indicating  if  the  user  wants 
to  save  the  value  of  this  function  for  tabulation. 

12.  DATANAME  -  if  it  is  decided  to  save  the  function  then 
what  type  function  is  it; 

F  -  the  function  itself, 

AI(x)  -  one  of  the  eight  analog  to  digital, 
EA(x)  -  one  of  the  eight  digital  to  anal:g. 

It  should  be  noted  that  functions  are  of  three  types.  If 
the  computer  finds  a  location  for  the  function  or  it  is  not 
hardware  dependent  then  the  function  is  described  as  "F  .  If 
the  function  is  one  of  the  other  two  types  "ai(x)  or  '  lA  .x)" 


then  they  represent  physical  fixed  registers  within  the 
computer.  If  the  litter  two  types  are  used  as  arguments  In  an 
equation  they  are  not  listed  as  inputs  under  NR; INPUTS  for 
the  function.  Since  they  are  physical  registers  they  are 
called  by  their  name  (ie.  A£(5))  and  not  by  lAP  or  IA£. 

D.  LOGICAL  3TRUCT0B2 

The  logical  structure  discussed  will  be  a  general  outline 
of  program  flow.  A  complete  discussion  is  covered  in  the 
example  problem.  Initially  the  control  designer  enters  the 
program  by  responding  "IBASE  I^ENU"  to  the  CP/M  operating 
system  prompt.  This  will  place  the  program  unaer  control  of 
MENU. CM!  which  initializes  the  program  and  enters  a  C»SE 
statement  that  provides  a  menu  selection  shown  in  Figure  15A. 
At  this  point  the  designer  usually  starts  a  new  problem  or 
recalls  an  old  problem  that  modifications  will  be  made  to. 
Notice  that  in  order  to  save  a  problem  the  user  must  select 
from  this  menu  option.  In  general  the  next  choice  of  the 
designer  is  to  enter  the  modification  phase  either  selecting 
to  change  only  the  initial  conditions  of  an  old  problem  -r 
selecting  to  make  major  changes.  If  major  or  new  changes  are 
chosen  the  program  transfers  control  to  CMENU.CMI)  where  the 
designer  is  allowed  to  change  the  function  database  by  menu 
selections  shown  in  Figure  16A.  The  choices  are  controlled  by 
the  CASE  statement  in  C^ENU.CME.  Most  selections  are  single 
command  entries  followed  by  a  ''RETURN'’  from  the  keyboard. 
Some  commands  will  operate  without  a  "RETURN'  especially  if 
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the  user  has  exceeded  a  field  in  one  of  the  database  records. 
The  fields  are  marked  on  the  display  with  a  at  each  end 
of  the  field.  On  some  terminals  the  arrow  keys  may  not  always 
work  for  positioning  the  cursor  in  a  field  desired.  In  this 
case  use  "CTIIL  E”  (control  key  pressed  at  the  same  time  as 
the  "e”  key)  to  move  the  cursor  up  a  field.  Use  "CTRL  X  to 
move  down  a  field  and  ’’CTRL  C"  to  move  down  a  complete 
record . 

Once  the  user  is  sure  the  database  contains  all  correct' 
functions  and  relationships  he  should  return  to  the  KAIN 
program"  if  he  wants  to  save  the  information  otherwise  any 
changes  will  be  lost  during  compiling.  The  program  could  be 
improved  by  adding  a  question  in  the  CONTROL. CMR  procedure 
asking  the  user  if  he  wants  to  save  the  information  before 
compiling.  It  takes  more  than  fifteen  minutes  to  compile  the 
program  so  a  few  extra  minutes  checking  for  correctness 
before  compiling  can  be  well  worth  the  effort. 

Once  the  option  to  compile  the  program  is  taken,  control 
passes  to  CONTROL. CMt.  This  is  the  train  control  procedure 
that  causes  the  automatic  programming  oy  calling  all 
necessary  procedures  required  to  write  the  PL/1  source  cole 
from  the  database.  A  complete  discussion  of  this  procedure  is 
contained  in  the  example  problem,  after  the  PL/1  source  .ode 
is  written,  a  control  program  is  written  on  how  to  compile 
the  ?L/1  code.  Finally  the  database  transfers  control  to  the 
control  program  it  Just  wrote,  called  FORM. SUB,  by  itlng 
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the  database  and  submitting  FORM. SUB  to  the  CF/M  operating 
system.  At  this  point  the  PL/1  programs  are  compiled  and 
executed. 


VI.  EESGRIPTIQN  CF  PL^l  ALGORITHM 


A.  GENERAL  DISCUSSION  (Refer  to  Appendix  B) 

This  section  describes  the  PL/1  programs  and  hew  they 
interact  to  provide  simulation  of  a  computer  controller  for 
an  analog  plant.  There  are  five  basic  proceaures  consisting 
of  a  monitor  that  insures  all  functions  are  calculated  at 
least  once  during  each  sample  period  of  the  computer,  an 
analog  procedure  for  describing  analog  calculations.  a 
digital  to  analog  procedure  and  analog  to  digital  proceiure 
for  describing  conversions,  and  a  aigital  p-ocedure  for 
describing  computer  calculations.  There  is  only  one  monit.^r 
but  there  are  as  many  of  the  other  procedures  as  needed  to 
describe  the  entire  system  model  for  each  function.  A:-  an 
example  if  there  are  two  analog  to  digital  equations  in  the 
problem  then  there  would  be  two  analog  to  digital  procedures, 
one  for  each  equation.  The  example  problem  will  derronstratr 
the  function  relationships.  The  basic  procedures  requ’ire 
information  for  each  function  and  that  is  nrovided  by  the 
control  engineer  throuth  the  database  program  in  the  fo'-^  of 
Include  files.  Once  the  program  is  ready  and  compiled 
automatically  by  the  database  program  a  result  file  called 
ANSwEH.COM  is  produced.  The  problem  is  run  by  typing  "aN/jiER' 
to  the  CP/M  operating  system  prompt.  As  the  program  is 
running  intermediate  results  are  displayed  on  the  lisvlay 


device  for  trou^Jleshootlng  should  the  results  not  he  as 
expected.  After  the  program  runs  it  produces  a  file  DATA.FLT 
containing  floating  point  data.  This  file  can  he  translated 
into  a  text  file  by  typing  FEINT  FILE  where  FILE  is  a  name  of 
a  file  provided  by  the  user  that  he  wants  the  results  placed 
into . 

PL/l  was  chosen  as  the  operating  system  language  because 
of  several  features.  First  it  is  a  common  high  level  language 
that  is  transportable  to  many  of  the  present  microcomputer's. 
The  intent  of  this  project  was  to  '<eep  the  programs  as 
machine  independent  as  possible  and  using  ?L/1  source  code  is 
a  good  way  to  do  this.  PL/l  also  has  very  good  math 
capabilities.  Since  the  type  of  problems  being  solved  are 
engineering  in  nature  it  was  necessary  to  have  a  high  level 
math  language.  Additionally,  the  aoility  for  PL/'i  to 
communicate  with  the  computer  hardware  is  a  definite  asset  i^ 
the  programs  are  to  be  realized  ir.tn  hardware.  Several  f'ther 
languages  were  lookea  at  such  as  PA.5CAL  but  PASCAL  dees  not 
have  the  ability  to  talk  to  the  hardware.  The  C  language 
was  also  looked  at  but  it  lacks  the  math  capability.  As 
stated  in  Chapter  IV. i  hopefully,  ACA  will  cverceme  all  the 
difficulties  by  combining  all  the  needed  features  into  one 
language . 

PL/1-60  has  some  problems.  There  is  a  definite  nroolem 
with  conversion  factors  as  was  discussed  in  Chapter  III.  The 
main  nroblems  have  to  do  with  conversions  between  floating 
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point  and  other  forms.  The  original  approach  was  to  use  only 
floating  point  values  and  therefore  not  need  conversions  but 
it  was  later  discovered  that  several  functions  do  conversions 
regardless  of  starting  and  ending  types.  This  is  a  prcblem 
and  cannot  be  overlooked  by  the  designer  if  he  has  problems 
with  results.  If  equations  requiring  ccmplei  functions  are 
used  the  designer  will  need  to  'oe:ome  familiar  with  PL/1. 

£.  LEVEL  #1  (v'AJO.a  5U1HCOTINE3  AMD  PEODEtJaiS) 

■jQ-SIIQBaPLI  “  t^ain  control  loop  program 

a.  Causes  Iteration  of  ail  functions. 

b.  Insures  all  functions  are  calleu  in  proper  crnar. 

c.  Saves  data  for  later  display  and  analysis. 

2*  ~  basic  format  for  all  analog  functions 

a.  Calls  all  functions  used  as  inputs. 

b.  Calculates  and  returns  it's  own  present  or  oast 
values  to  all  functions  that  call  it. 

c.  Insures  synchronization  of  calculations  with  the 
system  clocks  through  event  counts. 

“  basic  format  for  all  analog  tc  dl.-ital 
functions  and  performs  laeatical  functions  to  those  shown  in 
2a,  2b,  and  2c  above. 

"  ^^5ic  format  for  all  digital  functions  and 
performs  identical  functions  to  those  shown  In  2  above. 

£4.r££i  "  basic  format  for  all  digital  to  analog 
functions  and  performs  identical  functions  to  those  above. 
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6.  ”  corrects  all  improperly  formed  PL/1 
"include”  files  by  changing  the  backslash  produced  by  the 
BLASE  II  urogram  into  semicolons. 

7.  PBINT^P^I  -  utility  program  showing  how  to  convert  the 
floating  point  file  DATA.FLT  produced  by  the  main  program 
into  a  different  printable  text  file  for  display  on  the 
system  display. 

C.  LEVEL  #2  (PROGRAM  EESCRIPTION) 

Since  the  majority  of  the  code  is  in  higher  level 
languages,  and  is  designed  in  modules,  it  is  already  in 
algorithmic  form.  It  is  assumed  that  the  individual  needing 
to  understand  the  code  will  be  familiar  with  references  1 
through  4;  therefore,  an  extensive  discussion  of  the  computer 
code  will  not  be  provided.  The  example  problem  will  cover  all 
the  features  of  the  program. 

The  basic  structure  of  the  program  uses  a  monitor  that 
sequentially  calls  all  functions  within  the  computer  once 
during  each  digital  clock.  The  system  clock  is  used  as  an 
event  counter  to  provide  synchronization  for  each  function. 
The  functions  are  developed  by  the  control  engineer  through 
the  database  program  and  each  functional  relationship 
indicating  which  functions  call  each  other  must  be  uronerly 
provided.  This  Is  taken  care  of  when  the  designer  states 
which  functions  are  Inputs  to  each  other.  A  function  can  have 
its  output  connected  to  any  number  of  functions  but  a 
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functloa  can  only  have  five  inputs  not  counting  analog  to 
digital,  digital  to  analog,  and  it's  own  delayed  output.  The 
function  definition  allows  for  very  ccmplex  multiple 
input/output  relationships  within  any  single  equation.  A 
method  for  correcting  for  computation  and  other  delays  is 
provided  in  the  correction  factor  CF  in  the  monitor.  This  is 
a  factor  provided  for  the  erample  problem;  however,  it  is 
felt  that  this  equation  will  work  for  any  control  system  of 
the  type  requiring  a  steady  state  error  signal  of  zero.  The 
CF  is  provided  but  can  be  eliminated  simply  by  not 
referencing  it  in  any  of  the  equations.  The  programs  provided 
are  for  simulation  only  but  should  be  easily  realize-^  into 
control  hardware  since  the  design  equations  will  not  change. 


E.  LOGICAL  STRUCTURE 

The  program  contains  a  number  of  procedures  that  are 
basic  to  all  problems.  These  were  discussed  above. 
Additionally  each  problem  has  unique  features  that  are 
provided  by  the  control  designer  through  the  database 
programs.  These  features  are  coded  into  PL/1  "include”  files 
and  read  in  at  compile  time.  Once  the  programs  are  coded  all 
information  as  to  which  functions  are  connectec  together  will 
be  provided.  The  monitor  is  a  DO  #aiLE  LOOP  that  calls  all 
analog  to  digital  and  digital  functions  once  each  count  of 
the  clock.  At  each  count  of  the  digital  clock  [  DCLOCE  i  any 
called  function  will  perform  it's  own  equation  evaluation  if 
it  has  not  already  been  calculated  during  this  clock  and 
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return  the  correct  value  to  the  invoking  function.  When  the 
monitor  calls  a  function  It  does  not  use  the  value  returned. 
Additionally  there  Is  a  DO  WHILE  LOOP  Inside  the  main  monitor 
loop  used  to  call  the  analog  functions  the  number  of  times 
specified  by  the  designer  to  occur  between  each  digital 
sample  period.  This  Is  the  analog  clock  (  ACLOCK  ).  Also 
during  the  analog  clock  when  the  ACLOCK  is  equal  to  the  DELAY 
specified  by  the  designer  the  digital  to  analog  function  Is 
invoked.  When  a  function  Is  Invoked,  the  first  thing  that  It 
does  is  check  to  see  if  It  has  any  inputs.  If  there  are  any 
inputs  then  the  input  functions  are  called  to  Insure  they  are 
the  correct  value  before  the  equation  is  calculated.  Even 
though  the  monitor  is  a  sequential  process  invoking  each 
function,  once  a  function  Is  called  It  will  Invoke  all  of  the 
input  functions  needed  for  equation  calculation.  Some 
functions  could  be  called  many  times  during  each  clock 
depending  on  how  many  other  functions  use  it's  output.  Once 
each  ECLOCK  the  monitor  outputs  values  that  have  been 
requested  by  the  designer  to  save  into  a  file. 

The  first  procedures  called  by  the  monitor  are  all 
digital  to  analog  functions  located  in  a  file  written  by  the 
database  program  called  "ad.txt".  The  functions  take  the  form 
dummy_ptr  *  FDNCTION_NAMS  where  FUNCT I0N_NAME  Is  the  name  of 
the  function  and  the  file  FUNCTION_N'AME  contains  the 
necessary  information  to  represent  the  function.  The  example 
problem  has  a  more  detailed  discussion.  The  first  thing  the 
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digital  to  analog  functions  do  is  call  the  analog  functions 
that  are  used  as  inputs.  If  these  analog  functions  have  other 
analog  Inputs  they  call  their  inputs  before  calculating  their 
output  results.  This  process  will  continue  until  the  analog 
functions  have  no  inputs  or  their  inputs  come  from  a  digital 
to  analog  register  that  will  contain  a  value  for  a  correct 
point  in  time.  Similar  processes  occur  for  files  "dd.txt", 
"da.txt’*,  and  "aa.txt**  representing  digital,  digital  to 
analog,  and  analog  functions  respectively.  The  crder  of 
calling  is  controlled  initially  by  the  monitor  and  is 
extremely  important  for  all  processes  to  obtain  the  correct 
results.  Direct  coding  of  the  processes  would  be  much  more 
efficient  but  it  must  be  remembered  that  tnese  processes  must 
work  no  matter  how  the  designer  initially  defined  the 
problem.  This  method  of  problem  solution  has  merits  even  for 
the  experienced  programmer.  Dven  though  the  code  is  longer 
than  direct  programming  this  type  of  problem  will  typically 
require  one  or  two  hundred  files  introducing  the  possibility 
for  syntax  and  other  programming  errors. 

The  four  main  functions  aa.pli,  ad.pli,  dd.pli,  and 
da.pli  all  have  the  same  general  structure.  Basically  the 
function  is  called  by  the  monitor  or  another  function,  the 
function  calls  any  inputs,  and  then  returns  the  correct  value 
of  the  equation  at  the  time  the  equation  was  invoked.  This 
arrangement  places  a  process  on  the  stack  through  the  PL 
calling  routine  until  all  inputs  are  available  eliminating 
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the  need  for  developing?  a  special  process  ranager.  Usln^r  this 
approach  makes  the  programs  easily  adapted  to  concurrent 
processors  and  languages  such  as  AEA.  The  conversion  routines 
ad.pli  and  da.pll  do  not  have  a  provision  for  using  their 
delayed  values  since  they  would  not  he  available  In  real 
hardware.  This  forces  a  true  representation  of  the  hardware 
and  makes  the  designer  provide  Intermediate  holding  functions 
if  delayed  values  are  needed.  Additionally  the  actual  code 
for  hardware  will  be  easier  to  implement  at  a  later  date.  The 
other  functions  dd.pll  and  aa.pll  have  provision  to  return  a 
previously  calculated  value  if  only  the  delayed  value  is 
desired.  This  allows  the  code  to  run  much  faster  and  provides 
more  powerful  equations  within  a  single  function.  All  four 
processes  allow  for  up  to  five  function  Inputs  and  the 
procedure  knows  how  to  call  those  inputs  by  using  the 
variable  ''nr_lnputs"  to  determine  to  which  label  in  the 
function  to  jump.  Basically  all  functions  the  designer 
Implements  are  written  into  one  of  the  above  four  procedures. 
As  an  example  if  the  designer  specifies  a  digital  function 
called  DIGIT,  the  database  will  write  the  Include  files  and 
then  compile  the  dd.pll  procedure  into  a  file  ir.REL.  Next 
the  compiled  file  ED.REL  will  be  copied  into  DIGIT. R£L  wnich 
will  now  contain  all  the  necessary  information  fo’'  the 
function  DIGIT.  Additionally  the  necessary  information  such 
as  the  name  DIGIT  will  be  Included  into  the  monitor  so  that 
it  knows  to  call  the  function  once  each  clock  period. 


Before  any  of  the  PL/1  programs  can  he  compiled,  all  the 
text  files  written  by  EBASE  II  that  contain  a  backslash  must 
be  converted  to  semicolons.  This  is  accomplished  by  the 
procedure  REFORM. PLI.  This  is  a  simple  pro^tram  that  reads  in 
a  file,  translates  the  character,  and  outputs  to  a  new  file. 
The  files  reading  from  and  too  are  controlled  by  the  submit 
files  FORM. SUB  or  SFORM.3UB  written  by  the  EBASE  II  program. 

The  programs  compile  into  a  program  called  4\S'i(En.C0M 
under  the  control  of  the  database  prOf^,ram.  A\SkER.COM  is  the 
file  needed  to  run  the  problem.  The  designer  only  needs  to 
place  this  file  on  a  C?/M  based  system  and  type  ANSlEn  in 
response  to  the  operating  system  prompt.  This  will  run  the 
program  and  sand  the  intermediate  results  to  the  display  and 
send  the  floating  point  data  into  a  file  called  EATA.FLT. 

Once  the  program  has  been  run  it  may  be  desired  to  ulace 
the  floating  point  data  into  a  different  form  for  display  and 
analysis.  This  project  can  be  greatly  Improved  oy  providing 
the  results  on  a  graphical  display  device.  Both  tar.ular  and 
graphical  results  are  supplied  at  the  end  of  this  paper  and 
the  reader  can  determine  for  himself  which  is  the  best 
method.  One  utility  program  provlaed  Is  called  PRINT  PLI. 
This  program  reads  in  the  floating  point  data  by  using  the 
information  contained  In  the  file  DATAINFC.TXT  as  to  how  the 
data  was  stored  and  then  sends  the  information  into  a  text 
file  specified  by  the  user.  This  text  file  can  then  be 
printed  if  desired. 
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VII.  EXAMPLE 


A.  INTRODaCTION 

This  part  of  the  paper  will  make  reference  to  all 
aspects  of  the  project  and  provides  the  easiest  method  for 
understanding  the  entire  paper.  It  should  be  noted  that  the 
descriptions  as  to  how  the  data  is  entered  and  displayed  is 
somewhat  misrepresented  since  the  entire  project  was  designed 
for  an  interactive  environment.  The  actual  data  entry  at 
runtime  is  much  more  formal  with  data  fields  that  allow  only 
specific  entry  in  specific  areas  of  very  specific  type.  This 
is  not  a  constraint  but  an  aid.  is  an  example  if  the  user  is 
to  replay  if  he  wants  to  save  the  function  results  he  can 
only  answer  with  "T ,  ?,  Y,  or  N*’  for  "true,  false,  yes,  or 
no".  The  answers  can  be  in  upper  or  lower  case  and  i^  an 
improper  input  is  made  such  as  ”7"  the  computer  will  BEEP  and 
wait  for  a  correct  response.  Additionally  some  information  is 
supplied  in  reverse  video  or  highlighted.  In  the  event  that  a 
color  system  is  used  such  as  the  INTERCOLOR  the  information 
will  be  in  different  colors.  The  only  way  to  fully  appreciate 
the  aids  for  the  designer  is  to  run  the  program.  In  order  to 
provide  a  realistic  test  for  this  program  the  sample  problem 
will  include  a  plant  that  is  at  the  limit  of  stability  and  is 
initially  designed  to  be  controlled  by  a  bang  bang 
controller.  The  process  will  Include  a  control  loop  that 
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feeds  back  on  itself  which  provides  a  definite  headache  to 
computer  pro#;ratnirers  because  of  the  synchronization  problems. 


b.  PKOiLKM  SKTUP 
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Assuming  that  the  clamp  has  a 
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Figure  4A  ' 

Computer  and  Plant  Model 
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very  short  conversion  delay  compared  to  all  other  delays  la 
the  digital  system,  the  output  of  the  clamp  is  for  all 
practical  purposes  the  same  as  the  input  to  the  clamp.  This 
has  been  demonstrated  by  many  authors  [Ref.  5]  .  The  next 
initial  assumption  will  be  to  Ignore  any  delays  due  to 
conversions  and  calculations  in  the  computer  (ie  set  I  =  0). 
This  will  allow  for  simple  design  equations  yet  the  gross 
error  in  reasoning  will  be  easily  corrected  for  later.  '  It 
should  be  noted  that  most  control  systems  Ignore  the  delays 
due  to  computer  computation  since  the  delays  are  usually 
insignificant  compared  to  the  sample  rate.  i 

By  letting  X2(t)  =  the  plant  output,  Xl(t)  =  d/dttX2ft']. 
and  letting  the  input  to  the  plant  be  E2(t)  by  ignoring  the 
effects  of  the  clamp  and  computation  delay,  a  set  cf  discrete 
difference  equations  can  be  developed  describing  the  plant  at 
the  sample  instances  T.  These  equations  have  been  developed 
in  the  past  by  many  authors  [Ref.  6]  and  are  shown  bel'-w  in 
discrete  form: 

Xl(k  +  1)  =  Xl(k)  +  T  ♦  E2(lc)  (1) 

X2(k-*"1)  =  T  *  XKk)  +  X2(3c)  -  [(T  *  T)  /  21  *  E2(i;  i2^ 

where; 

d/dt[Xl(t)]  evaluates  to  Xl(k+1) 
d/dtfX2(t)]  evaluates  to  X2(lcfl) 

Now  the  next  task  is  to  control  the  plant  described  by 
the  equations  in  1  and  2  above  according  to  the  design 
requirements  for  a  minimum  time  response  ripple  free 
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controller.  Since  this  is  a  second  order  system  it  will  take 
two  time  periods  to  obtain  control.  Assume  the  Initial  state 


vector  equals  zero  (ie.  X1(0)  =  X2(0)  =  0).  Since  we  want  the 
second  time  state  with  Xl(2)  =  0  or  no  ripple  and  X2(2)  =  1, 

the  input  unit  step,  this  leads  to  the  set  of  equations  below 
after  some  manipulation: 

Xl(2)  =  T  »»  E2(0)  +  T  «  22(1)  =  0  (2^ 

X2(2)  *  [3*(T  *  T)/2]*E2(0)  +  [(T  *  T)/2]*22(l)  =  1  '4) 

Solving  equations  3  and  4  above  simultaneously  leads  to* 
22(0)  =  1/(T  *  T)  (5) 

22(1)  =  -  1/(T  *  T)  (6) 


By  definition  of  the  Z  transform,  equations  5  and  t  and 
the  fact  that  the  desired  output  requires  E2(k)  =  0  for  k  >  1 
leads  to: 

-1 

22(2)  =  1/(T  ♦  T)  -  fl/(T  *  T)]z  <7) 

Now  to  find  the  Z  transform  of  Si(k)  we  start  with  the 
problem  definition  of  zero  initial  conditions  and  a  unit  step 
input . 

At  k  =  0 


21(0)  =  U(0)  -  X2(0) 


=  1 

At  k  *  1 

0 

=  1 

EKl)  =  U(  1)  - 

X2(l) 

*  1 

f(T  * 

T)/2l  *  22(0) 

»  1 

[(T  * 

T)/2]  *  1/vT  *  T) 

»  1 

0.5 

=  0.5 
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At  k  »  2 

El(2)  =  0(2)  -  X2(2) 

=  1  -  1  =0 

For  K  >  1  El(k)  =  0 

The  above  set  of  equations  leads  to  the  Z  transfor'n; 

-1 

El(z)  *  1  +  0.52  (8) 

Since  E(z)  =  S2(z)  /  EKz),  the  algorithm  for  the 
computer,  we  can  now  take  the  ratio  of  equations  7  and  8  and 
cross  multiply  then  transform  to  find  the  difference 
equations  needed. 

-1  -1 

[l  +  0.5z  ]*E2(z)  =  [l/(T  «  T)]  *  [l  -  z  J^'SKz)  (5) 

E2(k)  +  0.5«E2(k-l)  =  [1/(T  ♦  T)]*[El(k)  -  El{k-l)]  10) 

S2(k)  =  -  0.5*E2(k-l)  +  [1/{T  *  T ) ] * [El ( k) -El ( k-1 ) ]  .11) 

Equation  11  is  the  one  used  for  the  computer  simulation 
ana  shown  in  Figure  4B  along  with  the  rest  of  the  model. 
Since  it  is  not  allowed  to  use  subscripted  or  arrays  for 
variable  names  in  this  program  implementation  the  following 
changes  have  been  made  in  Figure  4B  for  variables  in 
equations  1,  2,  and  11; 

SK  *  El(k) 

El(k-i)  referenced  as  the  delayed  value  of  EK  or  lAL  in 
the  function  E2K. 

E2S  *  S2(k) 

S2(k-l)  referenced  as  the  delayed  value  of  i2K  or  ”o  in 
the  function  E2K. 

XIN  *  Xl(k)  or  Xl(n)  if  n  different  from  k 


1 


DA(1) 


Xl(lc-l)  referenced  as  the  delayed  value  of  XIN  or  ”o"  In 
the  function  XIN  and  lAE  In  function  X2N1. 

X2N1  =  X2(k)  or  X2(n)  if  n  different  from  k 
X2(lc-1)  referenced  as  the  delayed  value  of  X2(k)  or  ”o” 
in  the  function  X2N1  and  lAI  in  function  AD2. 

UN  =  U(t) 

T1  =  [T  *  i]  /  j 

In  the  initial  design  the  sample  rate  of  the  digital 
system  T  and  the  sample  rate  of  the  analog  system  at  discrete 
points  T1  were  considered  the  same(  ie.  1  =  J).  The  only 
reason  for  providing  the  possibility  of  calculating  the 
analog  values  at  different  points  is  to  see  what  happens  to 
the  plant  between  sample  points  of  the  digital  system  at  T. 
In  most  problems  T  and  ”l”  will  be  normalized  to  unity  and 
"j"  will  represent  the  number  of  sample  points  calculated  by 
the  plant  for  each  unit  sample  period  T. 

Now  to  account  for  the  combined  delays  in  the  discrete 
system  due  to  analog  tc  digital,  digital  to  analog,  and 
computation  delays.  This  function  is  called  CF  or  the 
Correction  Factor.  The  CF  in  most  systems  is  a  ma.lor 
consideration  in  initial  equation  development  and  usually 
resembles  Kalman  Filters  or  predictors.  The  approach  taken 
for  this  thesis  was  that  the  average  designer  would  not  know 
or  could  not  estimate  these  delays  due  to  system  variables. 
This  means  the  simulations  would  not  work  as  expected  and 
needless  to  say  neither  would  the  realizations  of  the  control 
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algorithms  in  the  real  hardware.  The  attempt  was  to  force  the 

designer  to  consider  the  delays  in  his  design  and  possibly 

provide  a  good  predictor  initially  for  his  use.  It  was 

assumed  that  since  this  type  control  system  requires  an  error 

signal  that  is  zero  when  in  control  this  would  indicate  that 

the  signal  is  zero  mean.  Since  the  instantaneous  input  Is 

random  it  was  assumed  the  error  signal  was  gausslan  and 

provide  a  standard  exponential  predictor  to  the  error  signal 

-E 

to  correct  for  delays.  The  predictor  chosen  was  CF  =  e 
where  e  is  the  natural  log  and  E  is  the  combined  delays.  The 
C?  is  not  a  complete  guess  and  further  work  will  be  presented 
at  a  later  date.  The  results  of  the  CF  to  this  particular 
example  can  be  seen  in  the  plots  at  the  end  of  this  paper. 
Further  research  needs  to  be  done  to  see  if  this  correction 
works  for  truncation,  roundoff,  or  even  noise  in  the  input 
signal.  The  point  is  that  some  type  of  correction  needs  to  be 
added  and  this  one  works  very  well  in  this  case.  If  the  CF  is 
not  wanted  the  designer  only  needs  to  eliminate  it  from  any 
of  his  equations. 

C.  AUTOMATIC  PROGRAM  (Refer  to  Appendix  A  and  E) 

The  designer  will  Initially  enter  the  program  by  entering 
"EEAoE  menu”  and  a  RETURN  to  the  CP/M  operating  system 
prompt.  This  assumes  the  disk  containing  the  ERASE  programs 
is  in  drive  "a”  and  the  disk  containing  the  PL/1  programs  is 
in  drive  "E".  The  system  can  be  modified  to  run  on  a  single 
disk  as  was  done  In  the  original  system  but  the  disk  must 


contain  at  least  one  megabyte  and  should  have  directory  space 
for  at  least  one  hundred  fifty  files.  Ey  using  the  program 
setup  requirements  above  the  database  is  entered  by  menu 
selection  "l"  of  Figure  15A  to  start  a  new  program.  This  will 
initialize  the  system  and  enter  the  menu  selection  of  Figure 
16A .  Next  if  menu  selection  ”5"  is  entered  the  initial 
conditions  of  Figure  2  can  be  changed  to  the  ones  shown  in 
Figure  5A.  This  is  done  by  the  display  menu  CINIT  of  Figure 
17A.  Note  that  some  parameters  of  Figure  5A  are  effected  by 
returning  to  the  main  menu  and  saving  the  program  under  the 
name  STRT,  others  are  developed  as  part  of  the  program 
process  such  as  ACTION  which  changes  each  time  a  menu 
selection  is  made,  and  others  such  as  EATAl  are  developed  by 
the  program  itself  counting  the  number  of  items  that  have 
been  requested  to  save.  As  can  be  seen  from  Figure  5A, 
ALII^IT  was  entered  as  12  which  means  there  will  be  10  analog 
calculations  for  each  digital  calculation.  Since  the  initial 
design  used  the  same  sample  time  for  both  the  digital  and 
analog  then  Tl  must  be  adjusted  to  0.1  as  discussed 
previously  since  T  is  normalized  to  1.  After  the  initial 
conditions  ALIMIT,  DLIt'IT,  DELAY,  T,  and  Tl  are  entered,  the 
program  returns  back  to  CMENU  of  Figure  16A. 

Next  assume  menu  selection  "2”  is  selected  and  the 
functions  are  added.  The  structure  will  remain  basically  the 
same  with  only  slight  changes  shown  in  Figure  5B.  The 
database  for  the  problem  will  look  something  like  that  shown 
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in  Figure  6A  -  6E  belov..  Keep  In  mind  that  the  display  looks 
quite  different  due  to  color,  highlighting,  or  reverse  video. 
Also,  only  one  record  at  a  time  will  he  on  the  screen. 

The  remainder  of  the  menu  programs  will  not  be  covered 
since  they  are  self  explanatory.  It  should  be  noted  that  in 
order  to  save  the  problem  the  designer  is  working  on,  he  must 
first  return  to  the  main  menu  and  then  save  the  parameters. 
The  program  will  ask  for  the  name  of  the  file  the  user  wants 
to  save  the  program  in.  The  example  problem  was  saved  in  file 
3TRT.  If  there  are  questions  about  the  menus  look  at 
procedures  ^'EN^.C^'E  and  CMENU.CMi:  in  Appendix  A.  After  all 
functions  are  entered  CORRECTLY,  and  the  program  is  saved  if 
desired,  then  enter  “o”  for  menu  selection  in  Figure  16A. 
This  will  cause  the  program  to  enter  the  automatic 
programming  phase.  (Take  a  break  because  it  will  take  about 
fifteen  minutes.)  If  the  program  falls  to  compile  it  will 
probably  be  due  to  improper  naming  of  functions.  Function 
names  cannot  be  arrays.  If  the  results  are  not  as  expected, 
even  though  the  program  compiles,  then  the  equations  are 
probably  in  error  or  there  is  a  strange  conversion  in  FL/1. 
If  this  happens  then  unfortunately  references  2  through  4 
must  be  studied. 

Once  compiling  begins  control  passes  to  CONTROL. CME.  The 
first  function  called  is  ECLMCN.CME  which  writes  the  monitor 
declarations  ICLMON.TXT  for  MONITOR. PLI  and  files  AA.TTT, 
AI.TTT,  EE.TTT,  and  EA.TTT  containing  the  procedure  rails  for 
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RSCORL  00003 

PUNCNAME  :E2K  : 

FUNCTYPE  :Et; 

NR; INPUTS  :1: 

INPUT ;A  ;SK  : 

41P2r  :1: 

INPUTtB  ;TWO  : 

B1P2D  ;1: 

INPUT-.C  :  THREE  : 

P 1  P5)  p  •  1  • 

INPUTtD  :f6ur  : 

P 1 •  1  • 

INPUTiE  ‘.FIVE  ; 

ElP2r  ;l: 

FE  :  0.0000: 

FED  :  0  .0000 ; 

EQUATION  :F=(IAP-I4r  )/(m)-.5*0 
SAVEEATA  :T : 

DATA NAME  ;F  ; 


RECORD  00004 
FUNCNAME 
FUNCTYPE 
NR;INPUTS 
INPUT;  A 
A1P2B 
INPUT ;B 
BlP2r 
INPUT;C 
C1P2E 
INPUT  :E 
D1P2D 
INPUT  ;E 
ElP2r 
FE 
FEE 

EQUATION 

SAVEEATA 

EATANAME 


XIN  ; 

AA: 

0: 

ONE  : 

l: 

TWO  : 

1 : 

THREE  ; 

l: 

FOUR  : 

l: 

FIVE 
1 ; 

0  .0000 : 

0  .0000 : 
F=EA(l)’!'Tl-*-0 
T: 

F  : 


Figure  6B 

Database  For  Example  Problem 
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RECORD  00005 
FUNCNAME 
PUNCTTPE 
NR:INPDTS 
INPnT.-A 
A1P2D 
INPUT  :B 
B1P2E 
INPUT:C 
C1P2E 
INPUT:E 
D1P2D 
INPUT:E 
E1P2E 
FK 
FEE 

ECUATION 

3AVEDATA 

EATANAME 


X2NI  : 

AA: 

1: 

XIN  : 

2: 

TWO  : 

1 : 

THREE  : 

1: 

FOUR  ; 

l: 

FIVE  : 

1: 

0  .0000: 

0.0000: 

F=EA(l)*Tl*Tl/2+IAE*Tl+0 

T: 

F 


Figure  6C 

Eatabase  For  Example  Problem 
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RECORr  00007 

FUNCNAME 

:AD2  : 

FUNCTYPE 

:Ar: 

NR  .‘INPUTS 

:1: 

INPUT: A 

:E2N1  : 

A1P2E 

:2: 

INPUT.-B 

:TWO  : 

B1P2E 

:1 : 

INPUT  :C 

:THREE  : 

C1P2D 

:l: 

iNPUTtE 

:FOUR  : 

E1P2E 

:1; 

iNPUTtS 

tPIVE  : 

E1P2E 

•1  • 

*  X  * 

FK 

:  0  .0000  : 

FKD 

:  0.0000: 

EQUATION 

:Ar(2)=IAr 

SAVEEATA 

:T : 

DATANAME 

:AE(2) : 

RECORD  00009 

FUNCNAME 

EAl 

FUNCTYPE 

EA: 

NR:INPUTS 

1: 

INPUT: A 

12K 

A1P2E 

1: 

INPTJT:B 

TWO 

B1P2D 

1 : 

INPUT:C 

THREE  : 

C1P2E 

1: 

INPUT;E 

FOUR  : 

E1P2E 

1; 

INPUT-.Z 

FIVE  : 

E1P2E 

1: 

FK 

0 . 0000 : 

FEE 

0 .0000  : 

EQUATION 

DA(1)=IAP 

SAVEEATA 

T: 

EATANAME 

EAd): 

Figure  ec 

Database  For  Fxample  Problem 


the  rnonttor.  See  Figure  7.  Next  PROC.CMI  Is  Invoked  for  each 
of  the  four  types  of  functions.  This  procedure  writes  the 
declarations  for  each  function  called  ECLFDN.XXX  where  XXX 
stands  for  which  function  number  it  is  (ie.  A,  B,  C  etc.  and 
the  last  two  XX  indicate  the  function  type  (ie.  Kk,  AI,  EE, 
EA).  PROC.CME  also  writes  the  function  input  files  INPUT. XXX, 
function  equation  files  FUNC.XXX,  and  the  output  files 
OUT. XXX  that  are  used  to  indicate  which  element  in  the  data 
array  each  function  represents.  See  Figures  8A  through  8H. 
The  next  function  invoked  is  LATA.CML.  This  procedure  w.rites 
the  file  EATAINFC.TXT  that  tells  how  many  functions  are 
stored  in  the  floating  point  file  EATA.FLT  and  the  location 
of  each  element.  See  Figure  9.  Next  the  program  passes 
control  to  FORM.CME.  This  procedure  writes  the  program  that 
tells  the  computer  what  to  do  vith  all  the  files  that  have 
just  been  written.  Basically  it  produces  a  CP/M  SUBMIT  file 
that  tell  the  system  to  reform  or  correct  all  files  that  have 
a  backslash  into  files  containing  semicolons,  then  it  tells 
the  system  to  compile  each  function  and  place  it  under  it's 
own  unique  function  name.  After  compiling  each  function  it 
links  them  all  together  into  one  program  called  ANSWER. COM. 
See  Figure  Id.  Finally  after  the  database  program  writes  the 
SUBMIT  program  it  passes  control  to  the  program  it  Just  wrote 
and  creates  the  simulation  program  ANSWER. COM  as  just 
mentioned . 
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/*  ECLMON.TXT  */ 


UN 

XIN 

X2N1 

An 

Ar2 

DAI 

EX 

E2K 

T  floa 
T1  flo 
DLIMIT 
ALIf-IT 
delay 
1  dat 
2 


ezt 
ext 
ezt 
ext 
ext 
ext 
ext 
ext 
t  stat 
at  sta 
fixed 
float 
float 
ezt , 
data( 


entry{ fixed (7) ) 
entry( fixed (7) ) 
entry(fixed(7) ) 
entry( fixed (7) ) 
entry(flxed{7) ) 
entry(fixed(7)) 
entry(flxed(7) ) 
entry( f ixed (7) ) 
ic  ext  init( 
tic  ext  init( 
static  init( 
static  init( 
static  ialt( 


returns 
returns 
returns 
returns 
returns 
returns 
returns 
returns 
1.0000) , 
0.1000) , 
10) , 

10) , 

4), 


(ptr) 

(ptr) 

(ptr) 

(ptr) 

(ptr) 

(ptr) 

(ptr) 

(ptr) 


9)  dec(15,6) 


/*  AA.TXT  */ 

dutnfny_ptr  =  UN 
dur>n’y_ptr  =  XIN 
dummy'ptr  =  X2N1 

/»  AI.TXT  ♦/ 

duinmy_ptr  *  ADI 
dummy'ptr  *  AD2 

/♦  EE. TXT  */ 

dummy_ptr  =  EX 
durniny_ptr  *  E2X 

/*  EA.TXT  */ 

dummy_ptr  =  DAI 


(code)\ 
(code)\ 
(code  )\ 


( code  )\ 
(cede  )\ 


(cede  )\ 
( code  )\ 


(code )  \ 


Figure  7 

[Monitor  Eeclarations  and  Functions 
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/*  DCLPUN.AAA  */ 

UN  :proc(code)  returns {ptr)\ 

del 

nr^inputs  fixed(7)  static  tnlt(0), 

ONE  ext  eatry(f lxed(7) )  returns (ptr ) , 

TWO  ext  entry{ f lxed(7) )  Teturas(ptr ) , 

THREE  ext  entry(flxed(7) )  returns (ptr ) , 
roUR  ext  eatry(f lxed(7) )  returns (ptr ) , 

FIVE  ext  entryCf lxed(7) )  returns (ptr ) , 
cod(5)  fixed(7)  static  initd ,1 ,1 ,1 .1 ) . 

F  float  static  lnit(  1.0000), 

0  float  static  init(  0.0000)\ 
del 

1  dat  ext, 

2  data(  e)  deeds, 6)\ 

ICLINP.TXT  -  THIS  IS  PART  OF  ALL  FUNCTIONS  ANE  WAS 
FIXED  FOR  THE  THESIS  CASE  OF  NUMERIC 
CONTROL  SYSTEMS.  »/ 

del 

1  lA  Uasedda  ptr), 

2  lAP  float; 

2  lAE  float, 

1  IB  baseddb  ptr), 

2  IBP  float; 

2  IBE  float, 
lie  baseddc  ptr ), 

2  ICP  float; 

2  ICE  float, 

1  ID  baseddd  ptr), 

2  lEP  float; 

2  lEE  float, 

1  IE  based  (ie  ptr ) , 

2  lEP  float; 

2  lEE  float; 

/*  INPUT. AAA  ■*/ 
inputd)  =  ONE  \ 
lnput(2)  =  TWO  \ 
input(3)  a  THREE  \ 
inputlA)  a  FOUR  \ 
ln?ut(5)  a  FIVE  \ 

/*  FUNC.AAA  */ 

F»1.0000 

/*  OUT. AAA  */ 

data(  1)  a  ftc(?  )\ 

put  skip  list(''data(  l)  a  ',F  )\ 


Fissure  5A 

Function  Declarations  and  Equations 
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/♦  eclfun.baa  */ 


XIN  :proc(code)  returns (ptr )\ 

del 

Qr_lnputs  fixed(7)  static  init(0), 

ONE  ett  entry(ftied(7) )  returns (ptr ) , 

TWO  ext  entry(f ixed(7) )  returns(ptr) , 

THREE  ext  eatry(f lxed(7) )  returns(ptr ) , 

FOUR  ext  entry{f ixed(7) )  returns (ptr ) , 

FIVE  ext  entry(f ixed(7) )  returns (ptr ) , 
cod(5)  fiied{7)  static  init (1 , 1 , 1 ,1 . 1) , 

F  float  static  init(  0.0000), 

0  float  static  init(  0.0000)\ 
del 

1  dat  ext , 

2  data(  8)  dec(15,6)\ 

/*  ECLINP.TXT  -  THIS  IS  PART  OF  ALL  FUNCTIONS  AN!  WAS 

FIXEE  FOR  THE  THESIS  CASE  OF  NUMERIC 
CONTROL  SYSTEMS.  */ 

del 

1  lA  l3ased( la  ptr ) , 

2  lAP  floatT 
2  lAE  float, 

1  IB  based(ib  ptr), 

2  IBP  floatT 
2  IBE  float, 

1  IC  btised( ic_ptr  ) , 

2  ICP  floatT 
2  ICE  float, 

1  IE  based( Id  ptr  ) , 

2  IDP  floatT 
2  lEE  float, 

1  IE  based(ie  ptr), 

2  ISP  floatT 
2  lEE  float; 

/*  INPUT. BAA  ♦/ 
input(l)  *  ONE  \ 
input (2)  *  TWO  \ 
input(3)  »  THREE  \ 
lnput(4)  *  FOUR  \ 
input(5)  =  FIVE  \ 

/*  FUNC.BAA  */ 

?»BA(l)«Tl+0 

/*  OUT. BAA  */ 

data(  2)  =  ftc(F  )\ 

put  skip  list('data(  2)  =  ',F  )\ 


Figure  8i 

Function  Declarations  and  Equations 
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/*  rCLFUN.CAA  ♦/ 


X2N1  :proc(code)  returns (ptr  )\ 
del 

nr_inputs  fiied{7)  static  inlt(l), 

XIN  ext  entry( f lxed(7) )  returns (ptr ) , 

TWO  ext  entry(f iied (7) )  returns (ptr ) , 

THREE  ext  entry( f lxed(7 ) )  returns(ptr ) , 

FOUR  ext  entry(f ixed(7) )  returns (ptr  ) , 

FIVE  ext  entry( f lxed{7) )  returns ( ptr ^ , 
cod(5)  fixed(7)  static  ini t (2 , 1 , 1 , 1 , 1 ) , 

F  float  static  inlt(  0.0000), 

0  float  static  init(  0.0000)\ 
del 

1  dat  ext, 

2  data(  3)  dec(15,6)\ 

/*  CCLINP.TXT  -  THIS  IS  PART  OF  ALL  FUNCTIONS  ANE  WAS 

FIXED  FOR  THE  THESIS  CASE  0?  MUMEEIC 
CONTROL  STSTS^"S  .  ♦/ 

del 

1  lA  based(ia  ptr), 

2  lAP  floatT 
2  lAE  float, 

1  IB  baseddb  ptr), 

2  IBP  floatT 
2  IBE  float, 

1  rc  baseddc  ptr), 

2  ICP  floatT 
2  ICE  float, 

1  ID  based (ld_ptr  ) , 

2  lEP  floatT 
2  lEE  float, 

1  IE  based (ie_ptr  ) , 

2  I EP  float, 

2  lEE  float; 

/*  INPUT. CAA  ’*'/ 
input (1)  =  XIN  \ 
input (2)  *  TWO  \ 
inuut(3)  =  THREE  \ 
input(4)  =  FOUR  \ 
input(5)  =  FIVE  \ 

/■*  FUNC.CAA  */ 

F=EA(l)’*'Tl*Tl/2+IAE*Tl+0 

/♦  OUT. CAA  ♦/ 

data(  3)  =  ftc(?  )\ 

put  Skip  list('data(  3)  =  '  ,7  )\ 

Figure  6C 

Function  Declarations  and  Equations 


/<«  DCLFUN.AAD  V 


ATI  ;proc(code)  returns (ptr)\ 
del 

ar_lnputs  fixed(7)  static  init(l), 

U\~  ext  entry( f ixed(7) )  returns (ptr ) , 

TWO  ext  entry(fixed(7) )  returns (ptr ) , 

THREE  ext  entry(f ixed(7) )  returns (ptr ) , 

FOUR  ext  entry( f ixed(7) )  returns (ptr ) , 

FIVE  ext  entry(f ixed(7 ) )  returns(ptr ) , 
cod(5)  fixed(7)  static  init U  ,1 , 1 ,1 ,1 ) , 

F  float  static  inlt(  0.0000), 

0  float  static  inlt(  0.0000)\ 
id 

1  dat  ext, 

2  data(  8)  dec(15,6)\ 

/«  DCLINP.TXT  -  THIS  IS  PART  OF  ALL  FUNCTIONS  AND  WAS 

FIXED  FOR  THE  THESIS  CASE  OF  NUMERIC 
CONTROL  SYSTEMS  . 

del 

1  lA  basedda  ptr), 

2  lAP  float, 

2  lAD  float, 

1  IB  based(lb  ptr), 

2  IBP  floatT 
2  IBE  float, 

1  IC  baseddc  ptr), 

2  ICP  floatT 
2  ICE  float, 

1  IE  baseddd  ptr), 

2  IDP  floatT 
2  lEE  float, 

1  IE  based(le  ptr), 

2  lEP  floatT 
2  lEE  float; 

/*  INPUT. AAE  ♦/ 
input(l)  *  UN  \ 

input(2)  =  TWO  \ 
input(3)  =  THREE  \ 
input(4)  =  FOUR  \ 
input(5)  =  FIVE  \ 

/«  FUNC.AAE  ♦/ 

AD(1  )*IAP 

/*  OUT. AAE  */ 

data(  4)  =  ftc(Ar(l))\ 

put  skip  llst('data(  4)  *  '',AI(1))\ 


Figure  6E 

Function  Declarations  and  Equations 
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/♦  rCLrUN.BAE  */ 

AE2  ;proc(code)  returns (ptr )\ 
id 

nr_inputs  fixed(7)  static  Inlt(l), 

X2N1  ext  eatry(f lxed(7) )  returns (ptr ) , 

T'«0  ext  entry(f  Ixed  (7) )  returns  (ptr  ) , 

THREE  ext  entry{ f lxed(7) )  returns (ptr ) , 

FOUR  ext  entry(f ixed(7) )  returns (ptr ) , 

FIVE  ext  entry( f ixed(7) )  returns (ptr) , 
cod(5)  fixed(7)  static  init (2 ,1  *  1 »1 . 1 ) , 

P  float  static  init(  0.0000), 

0  float  static  inlt(  0.0000)\ 
del 

1  dat  ext, 

2  data(  a)  deeds, 6)\ 

/♦  ECLINP.TXT  -  THIS  IS  PART  OF  ALL  FUNCTIONS  ANE  WAS 

FIXED  FOR  THE  THESIS  CASE  OF  NUMERIC 
CONTROL  SYSTEMS.  '* / 

del 

1  lA  basedda  ptr), 

2  lAP  floatT 
2  lAE  float, 

1  13  baseddb  ptr ), 

2  IBP  floatT 
2  lEE  float, 

1  IC  baseddc  ptr), 

2  ICP  floatT 
2  ICE  float , 

1  ID  baseddd  ptr), 

2  lEP  floatT 
2  lEE  float, 

1  IE  based (ie  ptr  ) , 

2  lEP  floatT 
2  ISE  float; 

INPUT. BAD  */ 
inputd)  =  X2N1  \ 
input(2)  =  TWO  \ 
inputO)  =  THREE  \ 
input{4)  =  FOUR  \ 
input(S)  *  FIVE  \ 

/•f  FUNC.BAD  ■*/ 

Ad2)=IAL 

/*  OUT.BAE  */ 

data{  5)  *  ftc(AD(2))\ 

put  skip  list('data(  5)  =  ',AI(2))\ 


Figure  SE 

Function  Eeclarations  and  Equations 
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/♦  ICLFUN.AIE  */ 


£K  :proc(code)  returns (ptr)\ 

del 

nr  Inputs  fixed(7)  static  init(e), 

ONE  ext  entry( fixed(7) )  returns (ptr ) , 

TWO  ext  entry( f lxed(7) )  returns{ptr  ) , 

THREE  ext  entry( f ixed (7) )  returns (ptr ) , 

FOUR  exi  entry( f ixed{7  ) )  returns (ptr  ) , 

FIVE  ext  entryUlxed(7n  returns  (ptr ) , 
coa(5)  flxed(7)  statie  init (1 ,1 , 1 , 1 , 1 ^ , 

F  float  static  init(  0.(2000), 

0  float  static  init{  0.0000)\ 
del 

1  dat  ext, 

2  data(  3)  dec(15,6)\ 

/*  DCLINP.TXT  -  THI3  IS  PART  OF  ALL  FUNCTIONS  AN!  WAS 

PIXEL  FOR  THE  THESIS  CASE  OF  NUMERIC 
CONTROL  SYSTEMS .  */ 

del 

1  lA  based( ia  ptr  ) , 

2  lAP  floatT 
2  lAE  float, 

1  IB  based(ib  ptr), 

2  IBP  floatT 
2  IBL  float, 

1  IC  based ( ic  ptr ) , 

2  ICP  float, 

2  ICE  float, 

1  IE  based(id  ptr), 

2  IDP  float, 

2  lEE  float, 

1  IE  basedde  ptr), 

2  lEP  float, 

2  lEE  float; 

/*  INPUT. AEE  */ 
input (1 )  =  ONE  \ 
input  (?)  =  TWO  \ 
input (3)  =  THREE  \ 
input(4)  *  FOUR  \ 
input(5)  *  FIVE  \ 

/♦  FUNC.AEE  */ 

F=(AD(1)“AE(2))’S‘CF 

/*  OUT. AEE  */ 

data(  6)  =  ftc(F  )\ 

put  skip  list('data(  6)  »  )\ 


Figure  SF 

Function  Eeclarations  ana  Equations 
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/*  rCLFUN.BIE  */ 


E2K  ;proc(code)  returns (ptr )\ 
del 

nr_inputs  fixed(7)  static  init(l), 

EK”  ext  eatry(f ixed(7) )  returns (ptr ) , 

TWO  ext  entry( flxed(7) )  returns (ptr ) , 

THiEE  ext  etttry(fixed(7) )  returns(ptr ) , 

?OUH  ext  entry(f ixed(7) )  returns (ptr ) , 

FIVE  ext  entry(f ixed(7) )  returns(ptr) , 

cod(5)  fixed(7)  static  init (1 ,1 , 1 , 1 , 1 ) , 

F  float  static  init(  0.0000), 

0  float  static  init(  0.0000K 
del 

1  dat  ext, 

2  data(  a)  dec(15,6)\ 

/>!«  LCLINP.TXT  -  THIS  IS  PAST  0?  ALL  FUNCTIONS  AND  WAS 

FIXEI  FOB  THE  THESIS  SASE  OF  NUMERIC 
CONTROL  SYSTEMS.  '* / 

del 

1  lA  based( ia  ptr) , 

2  lAP  floatT 
2  lAD  float, 

1  IB  I)ased(lb  ptr), 

2  IBP  floatT 
2  IBD  float, 

1  IC  based{ ic  ptr ) , 

2  ICP  floatT 
2  ICD  float, 

1  IE  based( id_ptr ) , 

2  lEP  float, 

2  IDD  float, 

1  IE  basedde  ptr), 

2  lEP  floatT 
2  lED  float; 

/*  INPUT. BEE  */ 
input(l)  =  ES  \ 

in?ut(2)  =  TWO  \ 
input(3)  =  THREE  \ 
inputd)  =  FOUR  \ 
input(5)  =  FIVE  \ 

/*  rJNC.BEE  ♦/ 

?=(IAP-IAD)/(T’^T)-.5^0 

/*  OUT. BEE  */ 

aata(  7)  *  ftc(F  )\ 

put  ski?  list('data(  7)  =  ',F  )\ 


Figure  5G 
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tCLFUN.AIA  */ 


lAl  :proc(code)  returns (ptr )\ 

id 

nr  Inputs  flied(?)  static  Inlt(l), 

E2K  ext  entry{flxed{7) )  returns ( ptr ) , 

TkO  ext  entry(f Ixed (7) )  returns (ptr ) , 

THREE  ext  entpy( f ixed{7) )  returns (ptr } , 

FOUR  ext  entryCf lxed(7) )  returns (ptr ) , 

FIVE  ext  entry{f lxed(7) )  returns (ptr ) , 
cod(5)  fixed{7)  static  init ( 1 , 1 , 1 , 1 . 1 ) . 

F  float  static  init(  0.0000), 

0  float  static  init(  0.0000}\ 
del 

1  dat  ext, 

2  data(  8)  dec(15,6)\ 

/*  ECLINP.TXT  -  THIS  IS  PART  OF  ALL  FUNCTIONS  ANL  WAS 

FIXEI  FOR  THE  THESIS  CASE  OF  NUMERIC 
CONTROL  SYSTEMS.  */ 

id 

1  lA  basedda  ptr) , 

2  lAP  floatT 
2  lAU  float , 

1  IB  based ( ib  ptr ) , 

2  IBP  floatT 
2  IBD  float, 

1  IC  baseddc  ptr), 

2  ICP  float”; 

2  rCD  float, 

1  II  baseddd  ptr ) , 

2  IIP  float"; 

2  IDD  float , 

1  IE  based( ie  ptr)  • 

2  lEP  floatT 
2  lED  float; 

/*  INPUT. AEA  «/ 
inputd)  =  E2X  \ 
input (2)  =  TWO  \ 
input(3)  *  THREE  \ 
input(4)  »  FOUR  \ 
input(5)  =  FIVE  \ 

/*  FUNC.AEA  */ 

LA(1)*IAP 

OUT.  AC  A  '-V 

idta(  S)  »  ftc(EA(l))\ 

put  skip  list('data(  8)  =  ',IA(1))\ 


Figure  6H 

Function  ledarations  and  Equations 
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»'4V 


tAT\INJO.TXT 


/*  FORM. SUB  */ 

XSUB 

REFORM  AA.TXT  AA.TTT 
REFORM  AB.TXT  AD.TTT 
REFORM  El. TXT  ri.TTT 
REFORM  EA.TXT  EA  .TTT 
PLI  MONITOR 


REFORM  INPUT. AAA  INPUT. TTT 
REFORM  ECLFUN.AAA  DCLFUN.TTT 
PIP  FUNC.TXT=FUNC.AAA 
REFORM  OUT. AAA  OUT. TTT 
PLI  AA 

PIP  UN .REL=AA.5EL 

REFORM  INPUT. BAA  INPUT. TTT 
REFORM  DCLFUN.BAA  DCLFUN.TTT 
PIP  F0NC.TXT=FUNC  .BAA 
REFORM  OUT. BAA  OUT. TTT 
PLI  AA 

PIP  X1N.HEL=AA.REL 

REFORM  INPUT. CAA  INPUT. TTT 
REFORM  DCLFUN.CAA  DCLFUN.TTT 
PIP  FUNG. TXT*FUNC. CAA 
REFORM  OUT. CAA  CUT. TTT 
PLI  AA 

PIP  X2N1.REL=AA.REL 

REFORM  INPUT. AAD  INPUT. TTT 
REFORM.  ECLFUN.AAI  ECLJUN.TTT 
PIP  FUNC.TXT=*FUNC.AAL 
REFORM  CUT.AAE  OUT. TTT 
PLI  AE 

PIP  ADI .aEL=AE.REL 

REFORM  INPUT. BAE  INPUT. TTT 
REFORM  ECLFUN.BAI  ECLFUN.TTT 

PIP  func.txt=func.bad 
REFORM  CUT. BAE  OUT. TTT 
PLI  AD 

PI?  AE2.REL*AL.aEL 


Figure  1£A 
FORM. SUB 
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REFORM  input. Atl  INPUT  .TTT 
REFORM  UCLFUN.ALI)  UCLFUN  .TTT 
PIP  FUNG. TXT*FUNC. All 
REFORM  OUT  .Air  OUT. TTT 
PLI  DC 

PIP  Ee:.REL«rC.5EL 

REFORM  INPUT. BED  INPUT  .TTT 
REFORM  ICLFUN.ili;  tCLFUN.TTT 
PIP  FUNG. TXT=FUNC. BEE 
REFORM  OUT. BID  COT. TTT 
PLI  EE 

PIP  E2K.REL*EE.REL 

REFORM  INPUT.AEA  INPUT. TTT 
REFORM  ECLFON.AEA  EELFUN.TTT 
PIP  FUNC.TXTsFUNC.AEA 
REFORM  OOT.AIA  OUT. TTT 
PLI  EA 

PIP  EAl.R5L=DA.aiL 


PlI  PRINT 
LINK  PRINT 


LINK  A.NS«SR*MONlTOF.  ,UN  ,11N  ,X2N1  ,AD1  ,ADc  ,EA1  ,  IK  ,KZK 


Figure  105 
FORM. SUB 


E.  TEST  RESULTS  ANE  PROGRAM  EXECUTION 

As  mentioned  previously,  once  the  izwgzau  ANSWER. COM  is 
developed,  all  that  is  necessary  is  to  enter  ANSiEP.  to  the 
CP/M  operating  system  promnt.  Figure  4B  and  9  should  be 
referenced  during  this  discussion.  Figure  11  is  a  partial 
printing  of  what  will  be  on  the  display  when  the  program  is 
executing.  The  printout  is  for  the  sample  problem  but  with 
the  EELAy  set  to  zero  and  the  ELIMIT  shortened.  The  first 
thing  to  note  is  that  the  functions  are  not  executed  in  ' rder 
and  some  are  executed  several  times  during  each  deck.  Time  £> 
actually  runs  from  ECLOCK  =  0  to  ACLCCS.  =  1.  Ey  associ^.ting 
the  data  element  number  with  the  order  of  the  function  names 
in  figure  9  the  user  can  determine  which  function  is  being 
calculated  or  called.  As  an  example  the  first  element 
returning  an  answer  in  Figure  llA  is  element  4  which  eqi.ates 
to  All  in  Figure  9.  The  next  elements  are  £,  6,  and  7  that 
equate  to  AD2,  SX ,  and  S2K  respectively.  An  important  note  is 
that  since  there  is  no  delay,  function  3  (lAl)  or  the  di,-ital 
to  analog  function  is  calculated  or  called  at  time  0.  This 
will  not  be  the  case  when  discussed  later  on.  The  only 
function  that  has  a  value  is  element  1  or  JN.  There  is  no 
further  changes  in  the  data  until  DCLCCX  =  1  which  is  the 
same  time  as  ACLOCK  *  0  for  the  second  time.  At  this  point 
the  analog  to  digital  conversion  has  taken  place  for  the 
first  time,  the  computer  does  internal  calculation,  digital 
to  analog  is  complete  since  there  is  no  delay,  and  the  analog 
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plant  starts  to  integrate  the  signal.  This  is  really  the 
start  of  the  sample  since  the  computer  is  operating  on  the 
initial  conditions  for  the  first  ten  counts  of  the  analog 
clock.  By  the  first  sample  period  after  the  initial  values 
(ICLOCK  *  Z)  the  controller  output  changes  tc  -1  and  the 
plant  output  is  at  the  half  way  value  of  the  input.  ?y  the 
second  sample  period  after  the  initial  conaitior.s  'ECLOvii 
3)  the  plant  output  is  equal  to  the  forcing  function  i^nut 
ana  the  controller  output  is  0.  It  should  ne  noted  that  in 
all  programs  the  zero  point  of  the  sample  pe''iod  always 
starts  at  ECLOCtC  =  0  so  that  the  initial  conditions  cao  be 
graphed.  The  above  example  is  a  theoretical  non  reali-.alle 
simulation  of  a  control  system  with  no  calculation  delays. 
The  results  obtained  are  identical  to  classical  texti’^'clc 
exampl es . 

Figure  12  is  the  identical  probl-m  discussed  >.ith  the 
delay  equal  to  zero  but  it  shows  what  the  display  file  looks 
like  after  the  utility  program  PRINT. PLI  is  run.  The  results 
in  Figure  12  are  obviously  more  readable  sine®  the  function 
name  and  final  value  are  not  presented  until  after  all 
calculations  are  complete. 

Figure  13  is  a  partial  printing  cf  the  results  wren  csirs 
the  full  example  presented  above  with  a  forty  percent  delay. 
This  means  there  is  four  ACLOCK  counts  after  the  rcLCCK 
before  the  digital  to  analog  values  are  available.  The 
important  points  to  note  are  tha t  the  digital  t f  analog 


values  (DAI  or  function  8)  do  no  chants  until  the  analog 

clock  (  ACLOCK  )  is  equal  to  4  during  each  digital  clock,  the 

« 

overshoot  is  only  about  2.5  percent,  and  it  takes  five  times 
as  long  to  reach  stability  or  the  input  signal.  These  results 
are  much  more  representative  of  a  real  world  problem  and  the 
problems  that  would  be  encountered  with  a  true  hardware 
implementation . 
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VIII.  CONCLUSIONS 


A  major  point  of  this  thesis  was  to  demonstrate  the 
feasibility  of  writing  a  computer  pro^jram  that  can  be  used  on 
a  microcomputer  by  a  novice  programmer  to  develop  cor.trol 
system  models,  simulation,  and  realization.  This  type  c’ 
system  was  not  possible  in  the  past  due  to  system  memory 
limitations  and  unacceptable  software  programs.  The  approach 
taken  in  this  project  was  to  use  commercially  available  high 
level  programming  languages  to  reduce  the  coding  time.  As  was 
discussed  in  this  thesis,  the  languages  are  sot  sufficient  in 
themselves  and  more  than  one  lancua-'e  was  needed  to  ccnnlete 
the  project.  It  is  apparent  that  a  language  that  has  tice 
features  of  a  good  relational  database  couple'l  with  e-ctrer-ely 
good  mathematical  functions  and  the  acility  to  coTmun i 'a te 
with  the  computer  hardware  is  a  necessity  for  this  type  'f 
system  to  be  a  useful  tool. 

The  importance  of  recognizing  teat  control  systems  can  he 
modeled  and  manipulated  in  a  database  cannot  oe  oviTfiookel. 
Ly  carrying  this  approach  to  the  limit  it  is  possible  to 
develop  a  large  database  of  control  functions  that  could  be 
recalled  and  used  like  a  library.  Once  t-e  database  is  formed 
it  could  be  manipulated  by  artificial  Intelllrence.  The 
computer  Itself  could  te  told  to  develop  a  coctrcl  system  to 
match  a  desired  set  cf  waveforms  usln-i  the  information  stored 
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or  at  least  find  a  test  fit  solution.  Missing  parameters 
could  te  searched  for,  cr  best  fit  solutions  cculd  p-.ssltly 
be  provided  as  easily  as  usinf^  the  the  present  database 
systems  to  find  phone  numbers  and  names. 

It  is  significant  that  source  code  was  used  for  the 
intermediate  language.  A  special  purpose  program  couln  have 
been  written  that  provided  direct  fun;ition  calls  to  specific 
Hardware;  however,  by  using  high  level  language  sour'e  'ode 
the  program  will  run  on  many  different  machines  Just  by  usin^; 
their  resident  compilers.  This  iias  demonstrated  by 
transporting  the  proerams  to  four  different  machines  with 
different  operating  environments.  ly  using  source  code  the 
problem  was  kept  at  a  high  level  of  abstraction.  The  urograms 
the  computer  needed  to  write  were  similar  to  English  which 
helped  lessen  the  confusion  while  developing  more  than  one 
hundred  files  and  programs. 

k  fact  that  became  very  apparent  early  in  the  urogram 
development  was  the  need  for  a  correction  factor  to  ta-'e  into 
account  the  delays  in  the  system  due  to  calculations  and 
conversions.  This  problem  is  aadressed  in  many  models  bu;  it 
is  seldom  considered  in  most  computer  solutions  sin-'e  the 
computation  time  is  usually  insignificant.  It  was  metermine-i 
that  solving  for  all  possible  delays  was  next  to  Impossible 
especially  when  implemented  since  the  hardware  itself  would 
be  variable.  A  correction  factor  is  provided  that  works  very 
well  with  the. example  problem.  Furtner  investigation  needs  to 


be  done  usin^  other  functions  to  test  the  results;  however, 
it  Is  felt  that  the  factor  will  work  with  any  control  system 
designed  to  have  a  zero  steady  state  error  signal.  The  basic 
approach  is  to  have  the  designer  develop  the  model  assuming 
there  are  no  delays  and  then  estimate  the  delays  to  be 
applied  to  the  correction  factor.  If  the  factor  does  not  work 
then  try  to  develop  one  from  other  experience.  If  this  fails 
then  use  benchside  engineering,  make  estimates  and  wat:h  what 
the  results  do. 

Appendix  C  contains  several  graphs  of  data  covering  the 
example  problem.  A  variety  of  delays  were  introduced  am  the 
results  with  and  without  the  correction  factor  are  provided. 
The  problem  was  initially  designed  assuming  no  delay  and  taen 
corrected  for  with  the  correction  factor.  It  should  be  noted 
that  in  all  cases,  even  with  delays  up  to  £9  per  rent, 
acceptable  results  were  produced  within  five  times  the 
initial  sample  period  design.  This  is  significant  since  the 
uncorrected  functions  had  severe  problems  even  with  the  'mall 
delays.  4/hen  using  the  correction  factor  the  longer  delays 
corrected  better  and  it  is  felt  that  this  is  due  *o 
outputting  a  control  signal  closer  to  the  time  the  input  was 
going  to  sample  the  result  of  the  control.  If  a  signal  is 
output  and  sampled  immediately,  the  probability  that  the 
response  to  that  signal  is  the  one  observed  is  higher  than  if 
a  long  delay  were  encountered  possibly  resulting  la  severe 
data  error.  Sven  though  the  system  was  designed  for  a  step 


input  there  is  evidence  that  the  correction  factor  works  very 
well  for  a  ramp  input.  This  is  also  shown  in  Appendix  C. 

As  a  final  note  the  power  of  a  graphical  dynamic  display 
for  the  designer  is  invaluable.  Reading  through  many  pages  of 
data  makes  it  very  difficult  to  see  where  the  significant 
changes  are  being  made.  A  real  time  plot  will  cut  design  time 
by  an  order  of  magnitude. 
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CINIT  DBF 
CMENU  CMt 
CONTROL  CME 
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Programs  Required  for  Database  Operation 
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*  MEND. CMC 
SET  TALK  OFF 

RESTORE  FROM  INIT  ♦  See  Figure  2 

USE  I  NIT  *  See  Figure  2 

SET  FORMAT  TO  MENU  ♦  See  Figure  15A  and  151 

EO  WHILE  T 
READ 
DO  CASE 

CASE  ACTION  =  'o' 

ACCEPT  'DO  YOU  REALLY  WANT  To  OUIT?  Y/N  '  TO  CN 
IF  CN  =  'Y' 

CUIT 

SNDIF 
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COPY  TO  PARAMS 
USE  PARAMS 
DO  CMENU 
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ACCEPT  'ENTER  NAME  OF  OLD  PROGRAM  FILE  '  TC  PARAY 

RESTORE  FROM  ^PARAM 

USE  &PARAM 
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USE  PARAMS 
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ACCEPT  'ENTER  NAME  OF  OLD  PROGRAM  FILS  '  T'J  PARA/ 

RESTORE  FROM  &PARAM 
USE  SPARAM 
COPY  TO  PARAMS 
USE  PARAMS 
EO  SHORT 
ENECASE 
ENEEO 


file  MENU 


ENTE*T  0  TO  QUIT 

ENTER  1  TO  START  A  NEW  PROBLEM 
ENTER  2  TO  RECALL  AN  OLL  PROBLEM 
ENTER  3  TO  SAVE  PRESENT  PARAMATERS 
ENTER  4  TO  DISPLAY  PRESENT  FUNCTIONS 
ENTER  5  TO  CHANGE  INITIAL  CONDITIONS  ONLY 
(  This  will  allow  the  program 
to  compile  much  faster  ) 

ENTER  ACTION  #ACTION 


Figure  15A 
MENU 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


MENU.FMT 

7.17  SAY  "enter  0  TO  QUIT" 

8.17  SAY  ENTER  1  TO  START  A  NEW  PROBLEM" 

5.17  SAY  "ENTER  2  TO  RECALL  AN  OLL  PROBLEM" 

10.17  SAY  "ENTER  3  TO  SAVE  PRESENT  PARAMATERS " 

11.17  SAY  "enter  4  TO  DISPLAY  PRESENT  FUNCTIONS" 

12.17  SAY  "enter  5  TO  CHANGE  INITIAL  CONDITIONS  ONLY 
13,28  SAY  "(  This  will  allow  the  program" 

14,30  SAY  'to  compile  much  faster  )" 

16,23  SAY  "ENTER  ACTION" 

16,36  GET  ACTION 


Figure  15B 
MENU  Format 
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*  CMEMU.CME 

SET  FORf^AT  TO  C^1E^JU  ♦  See  figure  16A  and  163 
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CASS  ACTION  *  '3' 
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RECALL  RECORD  iCN 
CASE  ACTION  =  's' 

SET  FORMAT  TO  CINIT 
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file  CMENU  *♦* 


ENTER  0  TO  COMPILE  A  PROGRAM 

ENTER  1  TO  EEIT  TRANSFER  FUNCTIONS 

ENTER  2  TO  ALE  A  TRANSFER  FUNCTION 

ENTER  3  TO  DELETE  A  TRANSFER  FUNCTION 

ENTER  4  TO  RECALL  A  DELETED  TRANSFER  FUNCTION 

ENTER  5  TO  CHANGE  INITIAL  CONDITIONS 

ENTER  6  TO  RETURN  TO  THE  MAIN  PROGRAM 


ENTER  ACTION  FACTION 


Figure  16A 
CMENU 


CMENU. FMT 

9  3,15  say  "enter  0  TO  COMPILE  A  PROGRAM" 

9  4,15  SAY  "enter  1  TO  EEIT  TRANSFER  FUNCTIONS" 

{?  5,15  SAY  "enter  2  TO  ADD  A  TRANSFER  FUNCTION" 

(?  6,15  SAY  "ENTER  3  TO  DELETE  A  TRANSFER  FUNCTION" 

7,15  SAY  "enter  4  TO  RECALL  A  DELETED  TRANSFER  FUNCTION" 
0  8,15  SAY  "enter  5  TO  CHANGE  INITIAL  CONDITIONS" 

0  9,15  SAY  "enter  6  TO  RETURN  TO  THE  MAIN  PROGRAM" 

0  12,22  SAY  "enter  ACTION" 

0  12,35  GET  ACTION 


Figure  16B 
CMENU  Format 
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/'AD-A126  847  INTERACTIVE  MICROCOMRUTER  CONTROL  SYSTEM  MODELING  AND 
REALIZATION  USING  A  DATABASE  FOR  AUTOMATIC  PROGRAMMING 
(U)  NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CA  R  F  dOHNSRUD 
UNCLASSIFIED  DEC  82  F/G  9/2 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAl  BURtW  Of  SUNDARDS- I96J-A 


**♦*♦  File  CINIT 


LLIMIT  #DLIMIT 

ALIMIT  #ALIMIT 

lELAT  #EELAI 
T  #T 

T1  #T1 


(This  is  the  number  of  times  the  computer 
will  do  calculations,  ie.  The  number  of 
sample  periods. ) 

(This  is  the  number  of  times  each  analog 
calculation  will  be  done  during  each  sample 
period.  Effects  how  accurate  the  analog  plant 
is  modeled . ) 

(This  is  the  number  of  analog  calculations  done 
before  the  digital  to  analog  is  availible.) 

(This  is  the  sample  rate  of  the  digital  system. 
Tou  must  provide  it  if  you  are  going  to  use  it 
in  your  equations.) 

(This  is  the  sample  rate  of  the  analog  system. 

It  is  usually  the  same  as  T  above  during  design; 
however,  for  simulation  it  usually  works  out  to 
be  T  divided  by  ALIMIT  since  you  usually  want  a 
more  accurate  plot  of  the  analog  plant  to  see 
what  happens  between  the  digital  sample  periods. 
You  must  provide  Tl  if  you  use  it  in  your  equa¬ 
tions.)  (This  is  effectively  double  rate  sampling. 


) 


Figure  17A 
CINIT 


*  CINIT, 
(?  2,3 

0  2,10 
e  2,24 
0  3,25 

0  4,25 

0  5,3 

0  5,10 

0  5,24 

0  6,25 

0  7,25 

0  8,25 

0  9,3 

0  9,10 

0  9,24 

0  10,25 
0  11,  4 
0  11,10 
0  11,24 
0  12,25 
0  13,25 
0  14,  4 
0  14,10 
0  14,24 
0  15,25 
0  16,25 
0  17,25 
0  18,25 
0  19,25 
0  20,25 
0  21,25 


FMT 

SAT  "DLIMIT" 

GET  {LIMIT 

SAY  ^(This  Is  the  aumher  of  times  the  computer" 

SAY  "will  do  calculations,  le.  The  number  of" 

SAY  "sample  periods.)" 

SAY  ALIMIT" 

GET  ALIMIT 

SAY  ],’(Thls  Is  the  number  of  times  each  analog" 

SAT  ^calculation  will  be  done  during  each  sample" 

SAY  "period.  Effects  how  accurate  the  analog  plant" 

SAY  "is  modeled.)" 

SAY  EELAT 
GET  DELAY 

SAY  "(This  is  the  number  of  analog  calculations  done" 
SAY  "before  the  digital  .o  analog  Is  availlble.)" 

SAY  "T" 

GET  T 

SAY  "(This  is  the  sample  rate  of  the  digital  system." 
SAY  "You  must  provide  lt„lf  you  are  going  to  use  It" 
SAY  "in. your  equations.)” 
say  Tl 
GET  Tl 

Say  "(This  is  the  sample  rate  of  the  analog  system." 

SAY  ^^It  is  usually  the  same  as  T  above  during  design;" 
SAY  "however,  for  simulation  It  usually  worlis  out  to" 
SAY  "be  T  divided  by  ALIMIT  since  you  usually  want, a" 
SAY  "more  accurate  plot  of  the  analog  plant  to  see" 

SAY  ]|wfcat  happens  between  the  digital  sample  periods." 
SAY  "You  must  provide  Tl  if  you  use  it  in  your  equa-" 
SAY  tlons.)(Thls  is  effectively  double  rate  sampling.) 


Figure  17B 


CINIT  Format 


*  CONTROL. CM! 

EO  ECLMON 

STORE  0  TO  INI 

LOCATE  FOR  FUNCTYPE  =  'AA' 

EO  PHOC 

STORE  0  TO  INE 

LOCATE  FOR  FUNCTYPE  *  'AE' 

EO  PROC 

STORE  0  TO  INE 

LOCATE  FOR  FUNCTYPE  =  'EE' 

EO  PROC 

STORE  0  TO  INE 

LOCATE  FOR  FUNCTYPE  =  'BA' 

EO  PROC 

EO  DATA 

EO  FORM 

QUIT  TO  'SUBMITI  FORM' 

SET  TALK  ON 
RETURN 


se 


'-’a 


*  DCLMON.CMD 

SET  ALTERNATE  TO  B:ECL^'ON  *  B:  not  needed  If  single  disk 

SET  ALTERNATE  ON 

LOCATE 

10  WHILE  .NOT.  EOF 

?  '  '♦FUNCNAME+'  ext  en try ( f ixed (7 ) )  returns(ptr ) , ' 
CONTINUE 
ENECO 

COUNT  FOR  SAVEDATA  TO  DATAl 

?  '  T  float  static  ext  init( '’+str(T,8,4)  +  ') , ' 

?  '  T1  float  static  ext  init( '+str(Tl ,8,4)+ ' ) , ' 

?  '  ELIMIT  fixed  static  init('+  str(DLIMIT,4)+') , ' 

7  '  ALIMIT  float  static  init('’+  str(ALIMIT  ,4)+' ) , ' 

?  '  delay  float  static  lnlt('+  str(EELAY ,4)+ ' 

?  '  1  dat  ext,' 

?  '  2  data('+str(BATAl,3)-*-')  dec(15,6)' 

SET  ALTERNATE  OFF 
STORE  'AA'  TO  T:NAME 
DO  MONECL 

STORE  'AE'  to  TtNAME 
EO  MONECL 

STORE  'EE'  TO  TjNAME 
LO  MCNDCL 

STORE  'EA'  TO  T;NAME 

EO  MONECL 

RETURN 


♦  MONECL. CME 

SET  ALTERNATE  TO  B;&T;NAMS  *  B:  not  needed  if  single  disk 
SET  alternate  ON 
LOCATE  FOR  FUNCTTPE  »  T:NAME 
IF  EOF  THEN 
? 

ELSE 

EO  WHILE  .NOT.  EOF 

7  'dummy  ptr  »  '♦fuacname+ '  ^code ) \ ' 

CONTINUE” 

ENEEO 

ENEIF 

SET  ALTERNATE  OFF 
RETURN 


'  n 
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*  PROC.CMI 
10  WHILE  .NOT.  EOF 
STORE  IND+1  TO  INI 

STORE  'ECLFUN.'-*-$(F0N,INE.1)+FUNCTTPE  TO  T:NAME 

SET  ALTERNATE  TO  B:SiT:NAME  «  E:  not  needed  if  single  disk 

SET  ALTERNATE  ON 

?  FUNCNAME-*>':proc (code)  returns(ptr)\' 

?  'del' 

?  'nr_inputs  flxed(7)  static  init ( '+N3:INP0TS  +  ' ) , ' 

?  INPUT:A+'  ext  entry(f lxed(7) )  returns(ptr) , ' 

?  INPUT:B+'  ext  entry(f ixed (7) )  returns(ptr) , ' 

?  INPUT:C+'  ext  entry (fixed (7) )  returns (ptr) , ' 

?  INPUT:!*'  ext  entry(f ixed(7) )  returns(ptr ) , ' 

7  INPUT:E*'  ext  entry(f ixed (7) )  returns (ptr ), ' 

?  'cod(5)  fixed(7)  static  Inlt  ( '  +  lllP2r  +  ' , '+B1P2E+; 

' . '♦C1P2E+ ' , '♦E1P2I+ ' , '♦E1P2E+  ' ) ,  ' 

STORE  STR(FK,8,4)  TO  TEMPOUT 
?  'F  float  static  tnit( '+TEMPOOT  +  ' ) , ' 

STORE  STR(FKE,8,4)  TO  TEMPOUT 
?  '0  float  static  lnlt( '+TSMPOUT+')\' 

?  'del' 

?  '  1  dat  ext , ' 

?  '  2  datd( '♦STR{dat4l ,3)+')  dec(15,6)\' 

SET  ALTERNATE  OFF 

STORE  'INPUT. '*$(FUN, INt, 1 )*FUNCTYPE  TO  T:NAME 
SET  ALTERNATE  TO  B:&T:NAM£  *  B:  not  needed  if  single  disk 
SET  ALTERNATE  ON 
?  'inputd)  *  '♦INP0T:A*'\' 

?  'lnput(2)  *  '+INPUT:B+'\' 

?  'input(3)  =  '+INPUT:C*'\' 

?  'input(4)  »  '+INPOT;E+'\' 

?  'input(5)  *  '+INPOT:S*'\' 

SET  ALTERNATE  OFF 

STORE  'FUNC.'+$(FUN,INE,1)+FUNCTYPE  TO  T:NAME 

SET  ALTERNATE  TO  B:ST:NAME  *  B:  not  needed  if  single  disk 

SET  ALTERNATE  ON 

?  EQUATION 

SET  ALTERNATE  OFF 

STORE  'OUT.'+^(FUN,INr,l)+FUNCTYPE  TO  TiNAME 
SET  alternate  TO  B:RT:i\AME  ■*  B:  not  needed  if  single  disk 
SET  ALTERNATE  ON 
IF  SAVELATA 

STORE  ELEMENTS+1  TO  ELEMENTS 

?  'data( '+STR(ELEMENTS,3)  +  ')  *  f tc ( '+IATANAME+ ') V  ' 

7  Iput  Skip  llst( 'aata(]*STR(ELEMENTS,3)+L)  » 

77  IATANAME+')\' 

ELSE 

7 

ENEIF 

SET  ALTERNATE  OFF 
CONTINUE 
ENIIO 
RETURN 
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*  DATA.CMD 

SST  ALTERNATE  TO  R:£ATAINEO  E:  not  needed  if  single  disk 
SET  ALTERNATE  ON 

?  'timel  fixed  static  init( '+STH(ILIMIT*ALIMIT,4)  +  ') 

?  nr  data  fixed  static  lnit( '♦STR(EATA1,3)  +  'K' 

COaNT"FOR  SAVEDATA  .AND.  (FUNCTTPE  *  'AA')  TO  IND 
?  'nr  aa  fixed  static  lnlt( '+STR(IND,2)+'), ' 

COUNT  FOR  SAVEEATA  .AND.  (FUNCTTPE  *  'AE')  TO  INE 
?  'nr  ad  fixed  static  lnit{ '♦STR(IND,2)+') , ' 

COUNT  FOR  SATEEATA  .ANE.  (FUNCTTPE  *  'EE')  TO  INE 
?  'nr_dd  fixed  static  lnit( '■♦•STHdNE  ,2 )+ ') , ' 

COUNT  FOR  SATEEATA  .AND.  (FUNCTTPE  =  'EA ' )  TO  INE 
?  'nr  da  fixed  static  init( '+STH(INE,2)  +  ') , ' 

?  'naiiie( '♦STR(DATA1 ,3)  +  ')  char(6)  static  lnit(' 

STORE  'AA'  TO  T:NAME 
EO  EATAl 

LOCATE  FOR  FUNCTTPE  «  'AA ' 

IF  .NOT.  EOF 
?? 

ENEIF 

STORE  'AE'  TO  T:NAME 
EO  EATAl 

LOCATE  FOR  FUNCTTPE  *  'AE ' 

IF  .NOT.  EOF 

??  ',' 

ENEIF 

STORE  'EE'  TO  T:NAME 
EO  EATAl 

LOCATE  FOR  FUNCTTPE  »  'EE ' 

1?  .NOT.  EOF 

??  '/ 

ENEIF 

STORE  'EA'  TO  T:NAME 
EO  EATAl 
??  ')' 

SET  ALTERNATE  OFF 
RETURN 


*  DATA1.CMD 

LOCATE  FOR  SATEEATA  .ANE.  (FUNCTTPE  =  T;NAMS) 
EO  WHILE  .NOT.  EOF 
?  [']+FUNCNAME  +  ['] 

CONTINUE 
IP  .NOT.  SOF 
??  ',' 

ENEIF 

ENEEO 

RETURN 


lei 


*  FOHM.CME 

SET  ALTERNATE  TO  FORM. SUB 
SET  ALTERNATE  ON 

7  'B:'  *  B:  not  needed  if  single  disk 

7  'XSUB' 

7  'REFORM  AA.TXT  AA.TTT' 

7  'REFORM  AE.TXT  AE.TTT' 

7  'REFORM  BD.TXT  DD.TTT' 

7  'REFORM  lA.TXT  EA.TTT' 

7  'PLI  MONITOB' 

STORE  'AA'  TO  T:NAME 
EO  ECLFUN 

STORE  'AE'  to  T:NAME 
DO  DCLFUN 

STORE  'EE'  TO  T;NAME 
EO  ECLFUN 

STORE  'DA'  TO  T:NAME 
EO  ECLFUN 
7  'PLI  PRINT' 

7  'LINE  PRINT' 

7  'LINK  ANSWER-MONITOR' 

LOCATE 

DO  WHILE  .NOT.  EOF 
7?  ','+TRIM(FUNCNAME) 

CONTINUE 

ENDDO 

7 

SET  ALTERNATE  OFF 
RETURN 


♦  ECLFUN. CME 

LOCATE  FOR  FONCTTPE  -  TtNAME 
STORE  0  TO  INE 
EO  WHILE  .NOT.  EOF 
STORE  IND+1  TO  IND 

STORE  $(FUN,INE,1)+FUNCTYPE  TO  EXTtNAME 
EO  REFORM 
CONTINUE 
ENEEO 
RETURN 


*  REFORM.CMD 

7  'REFORM  INPUT. '+EXT;NAME+'  INPUT.TTT' 

7  'REFORM  ECLFUN. '♦EXT :NAME+'  ECLFUN. TTT' 

7  'PIP  FUNC.TXT-FUNC. '♦EXT; NAME 
7  'REFORM  OUT. '♦EXT; NAME*  '  OUT. TTT' 

7  'PLI  '+FUNCTTPE 

7  'PIP  '+TRIM(PUNCNAMS)+'.RBL-'+FUNCTYPE+'.fi 
RETURN 
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'4 


*  SHORT.CMD 


SET  F0RH4T  TO  CINIT  ♦  See  figure  174  and  17B 
READ 

SET  T4LC  OFF 

SET  ALTERN4TE  TO  B:DCLMON  *B:  not  needed  if  using  a  single  disk 

SET  ALTERNATE  ON 

LOCATE 

EO  WHILE  .NOT.  EOF 

?  '  '♦FUNCNAME+'  ext  entry (fixed (7 ) )  returns( ptr ) , ' 

CONTINUE 

ENLCO 

COUNT  FOR  SAVEIATA  TC  EATAl 

?  '  T  float  static  ext  init( '♦str (T,a,4)  +  ') , ' 

?  '  T1  float  static  ext  init  ( '+st r  (T1 ,8,4)-*- ' 

?  '  ELIMIT  fixed  static  lnlt{'+  str(ELIMIT,4)'*-') , ' 

?  '  ALIMIT  float  static  initi'*  str(ALIMIT,4)+') , ' 

?  '  delay  float  static  lnit('+  str(EELAT,4)+ ' 

?  '  1  dat  ext , ' 

7  '  2  dataC '♦str(EATAl.3)  +  ')  dec(l5,6)' 

SET  ALTERNATE  OFF 
SET  ALTERNATE  TO  SFORM.SUB 
SET  ALTERNATE  ON 
?  'XSUB' 

?  'B: '  *  B:  not  needed  if  using  a  single  disk 

?  'PLI  MONITOR' 

?  'PLI  PRINT' 

?  'LINK  PRINT' 

?  'LINK  4NSWER*M0NIT0R' 

LOCATE 

EO  WHILE  .NOT.  EOF 
77  ','+TRI«(FUNCNAME) 

CONTINUE 

ENEEO 

7 

SET  ALTERNATE  OFF 
QUIT  TO  'SUBMITI  SFORM' 

RETURN 


ie3 


APPENIIX  B 


PL/1-80  PROGRAMS 


AA 

PLI 

DD 

PL  I 

PL  10 

OVL 

PIP 

GOM 

REFORM 

3YM 

IE 

TXT 

FONC 

AAA 

FUNG 

BAA 

FONG 

GAA 

FUNG 

AAI 

FUNG 

BAB 

FUNG 

All 

FUNG 

3DD 

FUNG 

AIA 

AE 

ipfr 

INPUT 

TTT 

XIN 

REL 

AE2 

REL 

£K 

REL 

PRINT 

RSL 

ANSWER 

COM 

AI 

PLI 

LINK 

COM 

PLIl 

OVL 

REFORM 

PLI 

EGLMON 

TXT 

EA 

TXT 

OUT 

AAA 

OUT 

BAA 

OUT 

GAA 

OUT 

AAE 

OUT 

BAD 

OUT 

AIE 

OUT 

BDD 

OUT 

AEA 

EE 

TTT 

DCLFUN 

TTT 

UN 

REL 

All 

REL 

EA 

REL 

PRINT 

COM 

ANSWER 

SIM 

EA 

PLI 

MONITOR 

PLI 

PLI2 

OVL 

REFORM 

REL 

AA 

TXT 

ECLFUN 

AAA 

ECLFUN 

BAA 

ECLFUN 

CAA 

ECLFUN 

AAE 

ECLFUN 

BAE 

ECLFUN 

ADD 

ECLFUN 

BEE 

DCLFUN 

ADA 

•EAT  A  INFO 

TXT 

EA 

TTT 

CUT 

fp 

X2N1 

REL 

ED 

R  £i  L 

FUNC 

TXT 

PRINT 

STM 

PRINT 

PLI 

ECLINP 

m  y  rn 

PLI 

COM 

PLILIB 

lAL 

REFORM 

CCM 

AE 

TXT' 

INPUT 

AAA 

INPUT 

BAA 

INPUT 

CAA 

INPUT 

AAD 

INPUT 

BAD 

INPUT 

ADD 

INPUT 

3IE 

INPUT 

ADA 

AA 

TTT 

MONITOR 

nEL 

AA 

REL 

AE 

REL 

E2S 

REL 

EAl 

hSL 

DATA 

FLT 

Figure  13 

Programs  Required  for  FL/1-60  Operation 
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!r  V. 


MONITOR:  proc  optlons(maln«stack(512) ) ; 
del 

ftc  entry  (float  binary)  returns  (chard? )  var); 

del 

TIME  fixed  ext  static  init(0), 

DCLOCE  fixed  ext  static  init(0), 

4CL0CC  fixed  ext  static  init(0)« 

Ar(8)  float  ext  static  init(0,0,0,0»0, 0,2,0), 

DA(8)  float  ext  static  init (0,0, 0,0, 0,0, 0,0), 
duimny_ptr  ptr, 
code  fixed(7), 

CF  float  ext,  /*  correction  factor  for  error  signal  '^/ 
datafile  file; 

del 

%include  'dclmon.txt'; 

« 

9 

a 

open  f ile(dataf lie )  record  sequential  output  title( 'data. fit' ) 
env(b(5l2)); 

CF  *  f tc(exp(char(-delay/ALIMIT ) ) ) ; 
do  while  (DCLOCK  <«  DLIMIT); 

PUT  SKIP  LIST(ECLOCK,'rCLOC£' );  /*>  SIIT  */ 
code  =  i; 

ACLOCK  «  01 
%include  'ad.ttt'; 

%include  'dd.ttt'; 
code  «  i; 

do  while  (ACLOCK  <=  ALIMIT  -  1); 

PUT  SKIP  LIST(ACL0CK, 'ACLOCK', TIME, 'TIME');  EIIT  *  / 
if  ACLOCK  «  delay  then  do; 

^include  'da.ttt'; 
end;  /*  if  V 
%include  'aa.ttt'; 
write  f ile(dataf lie)  from(dat); 

TIME  «  TIME  ♦  i; 

ACLOCK  >  ACLOCK  +  i; 
end;  /♦  do  ♦/ 

DCLOCK  >  LCLOCK  i; 
end;  /*  do  ♦/ 
close  file(dataf lie); 
end  monitor; 


/*  AA.PLI  ♦/ 

%liiclude  'dclfun.ttt': 
del 

ftc  entry  (float  binary)  returns (cbar{17)  var); 
del 

TIM£  fixed  ext* 

T  float  ext, 

T1  float  ext, 

lnput(5)  variable  entryCfixed (7) )  returns (ptr) , 
EA  fixed  static  iait(0), 

EAl  fixed  static  Init(e), 

DA(8)  float  ext, 

ACLOCK  fixed  ext, 
code  flxed(7), 
result  ptr, 
in(5)  label, 

(la_ptr,lb_ptr,lc  ptr,id_ptr,ie  ptr)  ptr; 
^include  'dcllnp.txt  ; 

^include  'input. ttt'i 
if  code  «  1  then  do; 
if  EA  <  TIME  then  do; 

if  (nr_iuputs>  0)  Si  (ar_inputs<  6) 
then  do; 
m(l)»la; 
m(2)«lb; 
m (3  )»lc; 
iii(4)»id; 
m(5)»l'e; 

;50to  rn(nr_luputs ) ; 
end;  /*  then  ♦/ 
else  goto  If; 

le:  ie^ptr  «  input  (5)  (cod(5)); 

Id:  id^ptr  »  input  (4)  (cod(4)); 

Ic:  Ic^p^r  *  input  (3)  (cod(3)); 

lb;  lb_ptr  »  input  (2)  (cod  (2)); 

la;  ia  ptr  «  input  (1)  (cod(l)); 

If;  0  «  f; 

EAl  «  time; 

%inclnde  'func.txt'; 

EA  «  time; 
end;  /*  if  EA  V 
end;  /♦  if  code  then  ♦/ 
else 

if  EAl  <  TIME  then  do; 

0  ■  f; 

EAl  -  time; 
end;  /*  EAl  */ 

Xinclude  'out.ttt'; 
result  •  addr(F); 
return  (result); 
end; 
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/*  AE.PLI  ♦/ 

^Include  'dclfun.ttt 


t 


ftc  entry  (float  binary)  returns(char(l7)  var); 
del 

£A  fixed  static  init(0), 

T  float  ext, 

T1  float  ext, 

input(5)  variable  entry(fiied{7) )  returns(ptr ) , 
CF  float  ext , 

Ar(8)  float  ext, 
code  flxed(7), 
result  ptr, 

(ii(5)  label, 

( la_ptr,ib_ptr,ic_ptr,ld_ptr,ie_ptr )  ?tr; 

^include  'dclinp.txt'; 

Xinclude  'input. ttt'; 

if  EA  >  0  then  do; 

if  (nr_inputs>  0)  &  (nr_inputs<  3) 
thett“do; 
m ( 1  )»la; 
m(2 )»lb; 
m(3)*lc; 
m(4)*ld; 
m(5)*le; 

goto  m(nr_iaputs); 
end;  /•**  then  */ 
else  goto  If; 

le;  ie_ptr  =  input  (5)  (cod(5)); 

Id:  id_ptr  =  input  (4)  (cod(4)); 

Ic:  ic_ptr  =  input  (3)  (cod(3)); 
lb;  ib^ptr  =  input  (2)  (cod(2)); 
la;  ia’ptr  »  input  (1)  (cod(l)); 

If: 

^include  'func.txt'i  /*  AI(X)  =  lA  ■*  ETC 

« 

f 

end;  I*  then  ♦/ 
result  =  addr(F); 

EA  a  i; 

^include  'out.ttt'; 

f 

return  (result); 
end; 


/*  DD.PLI  */ 

Xlnclude  'dclfun.ttt'; 
del 

ftc  entry  (float  binary)  returnstchard?)  var); 

del 

TIME  fixed  ext, 

T  float  ext, 

T1  float  ext, 

input(5)  variable  entry(fiied (7 ) )  returns (ptr) , 
£r  fixed  static  inlt(0), 

EDI  fixed  static  init(0K 
CF  float  ext, 

AC(a)  float  ext, 

DCLOCK  fixed  ext, 
code  fixed (7) , 
result  ptr, 
m(5)  label, 

( ia_ptr,lb_ptr ,ic  ptr,id_ptr ,ie_ptr)  ptr; 
^include  'dcllnp . txt^J 

%include  'input. ttt'; 
if  code  *  1  then  doi 

if  EE  <  ECLOCE  then  do; 

if  (ar_inputs>  0)  &  dr_inputs<  6) 
then  do! 
m(l)*la; 
m(2 )=lb; 
m(3-)®lc; 
m(4)®ld ; 
m(5)»le; 

goto  m(ar_laputs ); 
end;  /*  then  */ 
else  goto  If; 

le:  le_ptr  =  input  (5)  (cod(5)); 

Id;  id^ptr  *  input  (4)  (cod(4)); 

Ic:  ic_ptr  =  input  (3)  (codlS)); 
lb;  ib_ptr  =  input  (2)  (cod  (2)); 
la;  id  ptr  =  input  (1)  (cod{l)); 

If;  0  =  f; 

EEl  *  ecloce; 

^include  'func.txt'; 

9 

EE  *  DCLOCK; 
end;  /♦  if  ED  */ 
end;  if  code  then  */ 
else  do; 

if  EEl  <  DCLOCK  then  do; 

0  *  f; 

EEl  »  eclock; 
end;  /*  EDI  ♦/ 
end;  /*  else  ♦/ 
iinclude  'out. ttt'; 
result  *  addr(F); 
return  (result); 
end ; 
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/*  DA.PLI  •/  . 

^include  'dclfun . ttt 

del 

ftc  entry  (float  binary)  returns (cha r( 17 )  var) 

del 

TIMS  fixed  ext, 

T  float  ext, 

T1  float  ext, 

Input(S)  variable  entry! fixed > 7) )  returns(ptr' 
CF  float  ext, 

AC(8)  float  ext, 

CA(8)  float  ext, 
code  flxed(7), 
result  ptr, 
m(f)  label, 

!la,ptr,lb.ptr,lc_ptr,ld_ptr,le_ptr)  ptr; 

Alnclude  'dclinp.txt'; 

^Include  'Input. ttt'; 

If  TIME  id  then  do; 

If  (nr_lnput5>  »?)  &  (nr_lnputsv.  6) 
then~(to; 
m !  1 )  ■  I  a ; 
m  ( 2 )  •  1  b ; 
m (3  )*lo; 
m!4) "Id ; 

(f  ( 5  )«le; 

«oto  m( nr_lnputs  ) ; 
end;  /*  then  */ 
else  goto  If; 

le:  le_ptr  ■  Input  \5)  (1^; 

Id;  Id'ptr  ■  Input  !4)  U); 

Ic:  lc_ptr  ■  Input  \3)  (1); 
lb:  lb~ptr  -  Input  (2)  (1); 
la :  la'ptr  ■  Input  (1 )  ( 1 ) ; 

If ; 

tlnclude  'func.txt';  /♦  IA(X)  -  IB  -  KCT 

t 

end;  /♦  then  ♦/ 
result  ■  addr(F); 

^Include  'out .  tt  t ' ; 

f 

return  (result); 
end; 


REFORM:  proc  options (main) • 
del 

(input, output)  file, 
buff  chard )  varying; 

open  file(input)  stream  env(b(l024))  title( '$1  .$1 ') ; 

open  file(output)  stream  output  env(b(1024))  title( '$2  .$2') ; 

do  while  (d'b); 

read  file  (input)  into  (buff); 
buff  ■  translate(buff , ' ; ' ,  A' ) ; 
write  file  (output)  from  (buff); 
end; 

end  reform; 


PRINT:  proc  options (main) ; 
del 

data  dec (15,6) , 

( j ,k,m)  fixed  bin  (?) , 

( time ,dat, tempos tr  )• char( 128 )  var , 

(info, input, output)  file, 

^include  'datainfo.txt'; 

t 

open  file(input)  record  sequential  title( 'EAT4.FLT ') 
env ( b(5l2) ) ; 

open  flle(output)  stream  output  tltle( 'd.$l') 
env(b(5l2) ); 

do  k  »  0  to  timei; 
time  »  '  '; 

temp_str»'TIME='l lchar{k); 
do  m  *  1  to  length( temp_str ) ; 

if  substr(temp_str ,m,I)  “*  '  '  then 
time  =  time  I  { sub5tr(  temp  5tr,iT;,l); 

time  ■  '  J  m  '  1 ! time; 

write  file(output)  from(time); 
do  J«  1  to  nr_data; 

read  fllednput)  lato(data); 

dat  ■  '^m*l  '  I  I  name  ( ,) )  1 1 ^ Ichar  (da  ta  ); 

write  file(output)  from(dat); 
end; 

end;  /*  do  k  ♦/ 
end; 
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AMPLITUDE 


2T  4T  6T  8T 

TIME 

Figure  15C 

Terminated  Ramp  Input  With  No  Corrections 


Figure  19D 

Terminated  Ramp  Input  With  Corrections 
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90%  DELAY 
UNCORRECTED 


90%  DELAY 
CORRECTED 


Figure  ISE 

Response  to  Varidble  3tep  Inputs  *'itb  and  Without  Corrections 


iS'V".  *  4^. 


...  -S'  V 


2 


90%  DELAY 
UNCORRECTED 


Figure  1?F 

fiesponse  to  Variable  Step  Inputs  with  and  Vitbout  Corrections 
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